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ABSTRACT 


A  Digital-Optical  Hybrid  Correlator  with  the  capability  of  pro¬ 
cessing  pulse  doppler  and  wideband  radar  signals  has  been  developed  based  on  the 
correlation  properties  of  an  optical  phased  array.  This  hybrid  system  retains  the 
flexibility  of  digital  technology  for  signal  storage  and  registration,  and  the  speed  of 
coherent  optical  technology  for  performing  the  signal  correlation.  Two  phased - 
array  modulators  have  been  developed.  The  Membrane  Light  Modulator  is  a 
mirror  whose  surface  can  be  deflected  in  localized  areas  by  electrostatic  forces. 
The  multichannel  crystal  modulator  employs  the  electro-optic  properties  of  lithium 
tantalate  to  achieve  phase  modulation.  Both  modulators  have  been  used  to  perform 
optical  phased-array  beam  steering.  Pulse  doppler  and  wideband  radar  signals 
have  been  processed  by  the  hybrid  correlator.  The  system's  capability  to  perform 
two-dimensional  correlation  of  frequency  diversity  radar  waveforms  has  been 
demonstrated. 


rrm 


1473 


UNCLASSIFIED 

Security  dlBSti  flea  lion 


UNCLASSIFIED 

Security  CUsftificatioa 


TG  1193  A 
SEPTEMBER  1972 


Technical  Memorandum 

A  DIGITAL-OPTICAL 
HYBRID  CORRELATOR 


R.  A.  MEYER 
D.  G.  GRANT 
J.  L.  QUEEN 


THE  JOHNS  HOPKINS  UNIVERSITY  •  APPLIED  PHYSICS  LABORATORY 
862 1  Georgia  Avenue  »  Silver  Spring,  Maryland  «  20910 

Operating  under  Contract  NQOOf  7?J?-C440f  with  the  Department  of  the  New 

■  1 1 

III 


Approved  for  public  release;  distribution  unlimited. 


IJJ,  I.  U  Jl  ,  ll.PUPIIJP* 


*  J'PM  * 1 !  IAI.,*W"-»W  IP  w » <>•  »■  UUP wv?  <wv-..n5yji*u  *f * 


rMC  JOHNS  HOSKINS  UNIVCMITV 

APPLIED  PHYSICS  LABORATORY 
Sn.vts  s»mw  Kwnjun 


CONTENTS 


List  of  Illustrations  ....  vii 

1.  Introduction  .....  1 

2.  Description  of  the  Optical  Phase  Modulators  .  3 

Membrane  Light  Modulator  ...  3 

Multichannel  Lithium  Tantalate  Phase 

Modulator  .....  11 

3.  Beam-Steering  Properties  of  an  Optical 

Phased-Array  Antenna  ...  19 

4.  Electronics  .....  29 

5.  One -Dimensional  Optical  Beam  Steering  .  33 

MLM  Beam  Steering  ....  33 

Crystal  Beam  Steering  ....  43 

6.  Two-Dimensional  Optical  Beam  Steering  .  53 

Optical  Configuration  for  Two-Dimensional 

Beam  Steering  with  the  MLM  .  .  53 

Selection  of  the  Element  Locations  in  the 

Thinned  MLM  Arrays  ...  55 

Two-Dimensional  Beam -Steering  Results 

Using  the  MLM  ....  57 

Two-Dimensional  Beam  Steering  Using 

the  Crystal  Modulator  ...  64 

7.  Correlation  Properties  of  an  Optical  Phased 

Array  .....  63 

8.  One-Dimensional  Radar  Signal  Processing  .  71 

Pulse  Doppler  Waveforms  ...  71 

Linear- Step  FM  Waveforms  ...  74 

Real-Time  Search  Mode  processor  .  .  75 

9.  Two-Dimensional  Radar  Signal  Processing  .  77 


Preceding  page  blank 


au  Jnwfiu' 


!TOff^|wr;.m'jiaiVi.;^gBMji»i^<yii n *  ■  ^■'.»,j 


TMI  JOHN*  MOM  INC  UMIVOWrTT 

APPLIED  PHYSICS  LABORATORY 

•k.v«k  OntiN*.  MMtnAMo 


References  .... 

Bibliography  .... 

Acknowledgments 
Appendixes 

A.  Deflection  Characteristics  of  the  Membrane 

Light  Modulator 

B.  Hybrid  Correlator  Electronics 

C.  Electro -Optic  Properties  of  the  Lithium 

Tantalate  Modulator 

D.  Optical  Phased  Array  Transform  Patterns 

E.  The  Effect  of  a  Small  Magnification  Change 

on  the  Focus  Position  for  a  Lens  near 
Unity  Magnification 

F.  Selection  of  the  Element  Locations  for  the 

Thinned  MLM  Array 


TMC  JOHN*  NOMIM  UMIVOMITV 

APPLIED  PHYSICS  LABORATORY 

»—IMO  MMVLANO 


ILLUSTRATIONS 


1  Cross-Sectional  View  of  the  Membrane 

Light  Modulator  ....  4 

2  Element  Array  Pattern  of  the  Membrane 

Light  Modulator  ....  5 

3  Cross-Sectional  View  of  a  Deflected  MLM 

Element  .....  7 

4  MLM  Imaged  onto  Transmit  ?*lask  at  Unity 

Ma  gni  f  i  ca  t  ion  .  .  .  .  .  8 

5  New  Versions  of  the  Membrane  Light 

Modulator  .....  12 

6  Multichannel  Lithium  Tantalate  Phase 

Modulator  .....  13 

7  Transmit  Mask  for  46  Channel  Crystal 

Phase  Modulator  ....  17 

8  Linear  Array  of  M  Equally  Spaced 

Isotropic  Radiators  ....  20 

9  Transform  Pattern  for  a  Linear  Array  of  6 

Equally  Spaced  Isotropic  Radiators  (M=6)  .  22 

10  Two-Dimensional  Optical  Phased  Array  in 

the  Front  Focal  Plane  of  a  Lens  .  .  24 

11  Grating  Lobe  Structure  for  M-by-N  Two- 

Dimensional  Array  of  Isotropic  Radiators  .  26 

12  Block  Diagram  of  the  Hybrid  Correlator 

Electronics  .....  30 

13  Schematic  of  the  Optical  Configuration  for 

One-Dimensional  Signal  Processing  with 

the  MLM  .....  34 

14  MLM  Transform  Pattern  ...  36 


-  vii  - 


THK  JOHN*  HOPKINS  UNIVtMITT 

applied  ph /sics  LABORATORY 


15  Central  Grating  Lobe  of  the  MLM  Transform 

Pattern  .....  38 

16  MLM  One-Dimensional  Beam-Steering 

Results  .....  39 

17  The  100  Rayleigh  Beam  Positions  of  the  MLM  .  41 

18  Intensity  at  Bore  sight  of  the  MLM  Transform 

Pattern  as  a  Function  of  the  Average  Phase 

across  the  MLM  Elements  ...  42 

19  Measured  Impulse  Response  Function  for  the 

MLM  Imaging  Lens  ....  44 

20  Optical  Configuration  for  One -Dimensional 

Signal  Processing  Using  the  Multichannel 

Lithium  Tantalate  Phase  Modulator  .  .  45 

21  The  Use  of  Cylindrical  Optics  to  Compress 

the  Illumination  Beam  in  One  Direction  .  46 

22  Transmit  Mask  for  46  Channel  Crystal 

Phase  Modulator  ....  47 

23  One -Dimensional  Beam-Steering  ResultB  Using 

the  Multichannel  Crystal  Modulator  .  .  49 

24  Rayleigh  Beam  Positions  over  Entire 

Ambiguity  Interval  for  the  46  Channel 

Crystal  Phase  Modulator  ...  50 

25  Grating  Lobe  Structure  of  the  Crystal 

Modulator  Transform  Pattern  .  51 

26  Two-Dimensional  Optical  Beam  Steering 

Using  the  MLM  ....  54 

27  Transform  Pattern  for  the  Three  Pseudo - 

Random  Transmit  Masks  ...  58 

28  Two-Dimensional  Beam-Steering  Results 

Using  a  Thinned  MLM  Array  .  .  61 

29  Conversion  of  a  Linear  Array  to  a  Two- 

Dimensional  Array  ....  65 


-  viii  - 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

SU.VIM  S»*INQ.  MMIUm 


30  Additional  Optics  Required  to  Perform  Two- 

Dimensional  Beam  Steering  with  the 

Crystal  Modulator  ....  66 

31  Pulse  Doppler  Signal  Processing  with  the 

Hybrid  Correlator  ....  72 

32  Block  Diagram  of  the  Receiver  for  Window 

Mode  Signal  Processing  with  the  Hybrid 

Correlator  .....  73 

33  Block  Diagram  of  a  Real-Time  Receiver  for 

the  Hybrid  Correlator  ...  76 

34  Two- Dimensional  Radar  Signal  Processing 

Using  the  Hybrid  Correlator  ...  79 


TMt  JOHNS  HOSKINS  UNIVfNWTY 

APPLIED  PHYSICS  LABORATORY 

■)LVIR  IMttH*.  kUMTUHO 


1.  INTRODUCTION 


Many  radar  environments  require  the  use  of  large 
time -bandwidth  transmissions.  To  process  these  wave¬ 
forms  by  strictly  electronic  means  is  prohibitive  because 
of  excessive  cost  and  size.  The  Fourier  transforming 
properties  of  a  coherent  optical  system  and  the  success 
of  photographic-film  optical  processors  in  mapping  opera¬ 
tions  have  stimulated  the  search  for  optical  modulators 
that  are  suitable  for  the  large  time -bandwidth  radar 
signal-processing  requirements  of  surface  weapon  con¬ 
trol  systems  where  processing  delays  must  be  held  to  the 
order  of  milliseconds.  Optical  processors  based  on 
ultrasonic  modulators  have  been  demonstrated,  but  are 
restricted  because  of  limitations  in  integration  time  (set 
by  propagation  time  in  the  acoustic  media).  Using  the 
correlation  properties  of  an  optical  phased  array,  a  hy¬ 
brid  digital-optical  system  has  been  designed  which  em¬ 
ploys  digital  shift  registers  for  the  signal  storage  and 
registration,  and  coherent  optics  for  the  two-dimensional 
correlation.  Such  a  system  maintains  the  flexibility  of 
digital  technology  and  the  speed  of  optical  technology. 

Two  high-speed  phased  array  modulators  have 
been  developed  under  the  auspices  of  the  Applied  Physics 
Laboratory  for  use  with  the  Digital -Optical  Hybrid  Cor¬ 
relator.  The  Membrane  Light  Modulator  (developed  by 
the  Perkin-Elmer  Corporation,  Norwalk,  Connecticut) 
is  a  mirror  whose  surface  can  be  deflected  in  localized 
areas  by  electrostatic  forces.  The  multichannel  crystal 
modulator  (fabricated  by  the  Isomet  Corporation, 

Palisades  Park,  New  Jersey)  employs  the  electro-optic 
properties  of  lithium  tantalate  to  obtain  phase  modulation. 
Both  modulators  have  been  used  successfully  in  a  one- 
dimensional  beam  steering  mode.  One-dimensional  cor¬ 
relation  of  pulse  doppler  and  linear -step  FM  waveforms 
has  been  successfully  achieved  using  the  hybrid  correlator 
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in  an  off-line  window  mode.  Two-dimensional  beam  steer¬ 
ing  with  the  Membrane  Light  Modulator  has  been  success¬ 
fully  performed  demonstrating  the  capability  of  performing 
two-dimensional  correlation  with  the  hybrid  correlator. 

While  the  system  has  been  operated  in  an  off-line 
mode,  it  can  be  directly  tied  into  a  suitable  radar  to  be 
used  in  the  real-time  window  mode.  The  two  electro¬ 
optic  modulators  and  current  digital  input  and  memory 
devices  are  suitable  for  use  in  on-line  processing  of 
search  mode  wideband  radar  signals  with  no  delay  other 
than  that  required  to  accumulate  the  echoes.  However, 
a  suitable  means  is  not  currently  available  to  output  the 
answers  form  such  a  two-dimensional  signal  processor. 
Suitable  devices  do  exist  for  one-dimensional  signal  pro¬ 
cessing  (i,  e. ,  pulse-doppler,  linear- step  FM). 
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2.  DESCRIPTION  OF  THE  OPTICAL  PHASE 
MODULATORS 


The  signal-processing  capabilities  of  an  optical 
phased  array  have  stimulated  the  search  for  high-speed, 
two-dimensional  optical  light  modulators.  Two  phased- 
array  modulators  have  been  designed  and  will  be  de¬ 
scribed  in  this  section. 


MEMBRANE  LIGHT  MODULATOR 

The  Membrane  Light  Modulator  (MLM)  is  a 
mirror  whose  surface  can  be  deflected  in  localized  areas 
by  electrostatic  forces.  A  diagram  of  the  physical  struc¬ 
ture  of  the  MLM  developed  for  the  Applied  Physics  Labo¬ 
ratory  by  the  Perkin-Elmer  Corporation  (Norwalk, 
Connecticut)  is  shown  in  Fig.  1.  The  complete  device 
is  fabricated  on  one  face  of  a  6 -mm -thick  optical  flat. 

/  thin  chromium -gold  film  is  vacuum  deposited  onto  this 
glass  substrate  and  then  etched  into  100  stripes  40  4m 
wide  and  on  50-4m  centers  to  form  the  electrode  struc¬ 
ture.  This  area  is  covered  with  a  thin  protective  layer 
of  silicon  monoxide  (SiO)  that  protects  against  electrical 
breakdown.  A  second  layer  of  SiO  about  1.  5  4m  thick  is 
applied  which  contains  10  000  individual  cylindrical  wells 
in  a  100-by-l00  array.  As  shown  in  Fig.  2, the  rows  of 
the  array  are  aligned  with  the  chromium-gold  stripes  of 
the  electrode  structure.  The  wells  are  37.  5  4m  in 
diameter  and  are  on  50-4m  centers  in  both  directions. 

The  next  layer  is  a  collodion  membrane  approximately 
0. 1  4m  thick.  The  membrane  is  coated  with  a  thin  layer 
of  antimony  which  serves  as  a  conducting,  reflecting 
surface. 

The  antimony  acts  as  a  common  electrical  ground 
for  all  of  the  elements.  When  a  voltage  is  applied  to  one 
of  the  electrodes  all  of  the  elements  in  the  row  deflect 
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ILLUMINATION 


Fig.  1  CROSS-SECTIONAL  VIEW  OF  THE  MEMBRANE  LIGHT  MODULATOR 


ELEMENT  ARRAY  PATTERN  OF  THE  MEMBRANE  LIGHT  MODULATOR 
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into  the  wells  due  to  the  electrostatic  force.  The  shape  of 
the  deflected  modulator  elements  is  approximately  para¬ 
bolic  (Fig.  3)  and  the  phase  of  the  light  reflected  off  an 
individual  element  can  be  described  by: 

(1) 


where  <pQ  is  the  peak  phase  reflected  off  the  center  of  a 
modulator  element,  r  is  the  polar  coordinate  originating 
from  the  center  of  the  element,  and  RQ  is  the  radius  of 
ithe  element. 

In  order  to  approach  flat  or  "piston -like"  phase 
modulation  (and  also  to  eliminate  the  light  reflected  off 
the  area  between  the  elements),  the  parabolic  elements 
are  masked  such  that  only  the  light  from  the  central  area 
of  the  modulator  elements  is  used.  This  is  accomplished 
by  imaging  the  modulator  (at  unity  magnification)  onto  an 
opaque  mask  that  consists  of  a  two-dimensional  array  of 
circular  apertures  located  on  the  same  centers  as  the 
modulator  elements  (Fig.  4). 

Assuming  perfect  imaging,  the  average  phase 
across  a  mask  aperture  can  be  determined  by  the  equation: 


where  is  the  average  phase  across  the  mask  apertures, 
R0  is  the  radius  of  the  modulator  element,  and  R  is  the 
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Fig.  3  CROSS-SECTIONAL  VIEW  OF  A  DEFLECTED  MLM  ELEMENT 
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Fig.  4  MLM  IMAGED  ONTO  TRANSMIT  MASK  AT  UNITY  MAGNIFICATION 
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radius  of  the  mask  apertures.  As  the  aperture  diameter 
is  decreased,  the  average  phase  (<0^  y)  approaches  the 
peak  phase  (<pQ)  and  therefore  the  light  transmitted  through 
the  mask  appears  more  like  a  piston  modulator.  However, 
the  total  power  decreases  as  the  aperture  diameter  de¬ 
creases.  Asa  compromise  between  these  two  factors, 
20-^m-diameter  apertures  are  used  (20  #m  is  approxi¬ 
mately  one-half  the  diameter  of  the  modulator  elements). 
The  average  phase  across  the  20-^tm  apertures  is: 


(3) 


The  peak  phase  reflected  off  a  modulator  element 
is  a  function  of  the  voltage  applied  to  the  electrodes.  For 
small  deflections,  the  peak  phase  is  approximately  pro¬ 
portional  to  the  square  of  the  voltage  applied: 


<p  «  kV 
o 


(4) 


where 


2  tr 
X 


cR 


8TD 


c  =  the  permittivity  of  the  air  in  the  well, 

R  =  the  radius  of  the  well, 
o 

D  =  the  depth  of  the  well, 

T  =  the  tension  of  the  membrane  (assumed 
to  be  constant  for  small  deflections), 
and 

X  s  the  wavelength  of  the  illuminating  beam. 


However,  for  larger  deflections,  the  peak  phase  is  a  func¬ 
tion  of  higher-order  powers  of  voltage.  A  more  detailed 
description  of  the  theoretical  relationship  of  phase  and 
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voltage  is  given  in  Appendix  A.  The  modulator  electronics 
have  been  fabricated  under  the  assumption  that  these  higher 
order  terms  are  negligible. 

When  overvoltages  are  applied  to  the  modulator 
elements,  they  "bottom  out. "  The  collodion  membrane 
adheres  to  the  thin  SIO  protective  layer  over  the  electrode 
structure.  The  maximum  stable  deflection  is  44%  of  the 
initial  separation  <D).  If  it  is  desired  to  utilize  a  25%  de¬ 
flection,  a  further  increase  of  only  10%  in  voltage  will 
cause  instability  (see  Appendix  A).  Tl.e  modulator  ele¬ 
ments  were  also  found  to  bottom  because  of  the  creepage 
of  the  membrane  when  moderate  voltages  were  applied  for 
long  periods  of  time.  The  electrical  protection  circuits 
described  in  Appendix  B  were  added  to  prevent  the  modu¬ 
lator  elements  from  bottoming. 


In  order  to  determine  the  time  response  character¬ 
istics  of  the  modulator  elements,  the  membrane  stretched 
over  the  cylindrical  wells  can  be  considered  as  a  drum 
head  having  natural  frequencies  of  vibration.  The  funda¬ 
mental  resonant  frequency  is  given  (Ref.  1)  by: 


f  = 


2.4 

2ffR 

o 


(5) 


where 

R  =  the  radius  of  the  wells, 
o 

T  =  the  tension  of  the  membrane,  and 
o  ~  the  mass  of  the  membrane  per  unit  area. 

In  a  vacuum,  the  MLM  surface  resonates  at  about  2.  5 

with  a  Q  of  about  30.  If  air  is  trapped  in  the  wells, 
the  system  appears  stiffer  because  of  air  compression 
and  the  resonant  frequency  is  raised  to  about  3  MHz. 
However,  the  vibration  is  more  strongly  damped. 


TM«  JOHNS  HOSKINS  UNIVOtSmr 

APPLIED  PHYSICS  LABORATORY 

■NLVm  SWWWK.  HkMUW 


Since  the  membrane  permits  air  to  diffuse  through 
it,  there  will  be  a  slow  relaxation  of  the  membrane  under 
a  steady  deflecting  force  as  the  air  escapes.  The  time 
required  for  the  membrane  to  reach  a  steady-state  posi¬ 
tion  is  in  the  order  of  25  to  100  ys.  Because  of  this,  the 
modulator  is  operated  in  a  vacuum  when  fast  time  re¬ 
sponse  is  required. 

Several  variations  of  the  MLM  have  been  designed 
by  Perkin-Elmer  on  internal  research  and  development 
(R&D)  programs.  In  the  self-masked  modulator  the  trans¬ 
mit  mask  is  an  integral  part  of  the  modulator  (Fig.  5a). 
The  illumination  is  in  the  opposite  direction  of  our  MLM 
and  circular  holes  have  been  etched  in  the  opaque  elec¬ 
trodes  to  provide  the  masking.  This  MLM  does  not  re¬ 
quire  the  use  of  an  imaging  lens  (and  therefore  avoids  all 
the  alignment  problems  associated  with  it).  In  the  Elec¬ 
tron  Beam  MLM,  the  membrane  of  the  self-masked  modu¬ 
lator  is  coated  with  small  circular  pads  —  one  for  each 
modulator  element  (Fig.  5b).  An  electron  gun  "writes” 
a  charge  density  on  the  pads.  Since  the  charge  at  each 
element  can  be  controlled  independently,  this  version  is 
an  inherently  two-dimensional  device  (Ref.  2). 


MULTICHANNEL  LITHIUM  TANTALATE  PHASE 
MODULATOR 

In  the  past,  the  high  voltages  required  for  2 ff 
phase  modulation  have  limited  the  use  of  electro-optic 
crystals.  This  section  describes  a  multichannel  lithium 
tantalate  phase  modulator  whose  speed  and  2rr  voltage  are 
compatible  with  high-speed  microelectronics. 

A  multichannel  lithium  tantalate  phase  modulator 
(Fig.  6)  has  been  fabricated  for  the  Applied  Physics  Labo 
ratory  by  the  Isomet  Corporation  (Palisades  Park,  New 
Jersey).  The  crystal  is  cut  and  polished  such  that  the 
surfaces  are  perpendicular  to  the  X^,  Xo,  X3  axes  (cor¬ 
responding  to  the  S,  b,  and  b  axes  of  the  crystal). 
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'Jhe  bottom  surface  of  the  crystal  (perpendicular  to  the 
X3  axis)  is  coated  with  a  thin  layer  of  chromium-gold 
which  acts  as  a  common  electrical  ground.  This  surface 
is  then  bonded  to  a  brass  block  whose  thermal  expansion 
coefficients  are  approximately  equal  to  those  of  lithium 
tantalate  so  that  strain  is  not  induced  in  the  crystal  due 
to  thermal  expansion.  The  brass  block  also  acts  as  a 
mechanical  clamp  to  inhibit  changes  in  the  crystal  refrac¬ 
tive  indices  due  to  the  piezoelectric  and  photoelastic  ef¬ 
fects  (Ref.  3).  A  thin  layer  of  chromium-gold  is  deposited 
on  the  top  surface  of  the  crystal  and  etched  into  46  parallel 
electrodes  0.  2-mm  wide  and  on  0.  5 -mm  centers.  These 

A 

electrodes  are  parallel  to  the  X2  axis  of  the  crystal. 

The  electro-optic  properties  of  single  domain 
lithium  tantalate  have  been  well  documented  (Ref.  4). 
Lithium  tantalate  is  uniaxial  and  remains  uniaxial  (with¬ 
out  a  rotation  of  the  principal  axes)  when  an  electric 
field  is  applied  parallel  to  the  X3  axiSA(Ref.  a6).  Coherent 
light  polarized  parallel  to  either  the  Xj  or  X3  axes  and 
propagating  along  the  X2  axis  will  remain  linearly  polar¬ 
ized  Aas  an  electric  field  is  applied  which  is  parallel  to 
the  X3  axis.  The  crystal  exhibits  birefringence  (both 
natural  and  induced)  and  therefore  light  which  is  not 
polarized  parallel  to  the  Xj  or  X3  axis  will  in  general 
become  elliptically  polarized. 

When  a  voltage  is  applied  to  one  of  the  electrodes 
of  the  phase  modulator  described  above,  the  electric 
field  under  the  central  area  of  the  electrode  is  approxi¬ 
mately  parallel  to  the  X3  axis.  Assuming  that  the  field 
is  uniform,  the  light  which  propagates  under  the  elec¬ 
trode  will  be  retarded  in  phase  (due  to  the  linear  electro¬ 
optic  effect)  by  an  amount: 


"ff  l2  n3  r33  V3 


(6) 
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A 

for  light  polarized  parallel  to  the  X3  axis,  and: 
-ff  l2  n?  r!3  V3 

iA<p3~ - u: - - 


for  light  polarized  parallel  to  the  Xj  axis,  where 

f2  =  the  length  of  the  crystal  in  the  X2  direc¬ 
tion  (15  mm), 

A 

=  the  thickness  of  the  crystal  in  the  X3 
direction  (0,  10  mm), 

n3  =  the  natural  index  of  refraction  for  light 
polarized  parallel  to  the  X3  axis  (2,  180), 

n^  =  the  natural  index  of  refraction  for  light 
polarized  parallel  to  the  axis  (2.  175), 

r  =  a  linear  electro-optic  coefficient 
(3.0  x  10-9  cm/V), 

r  -  a  linear  electro-optic  coefficient 
(6.  9  X  10-10  cm/V), 

X  =  the  wavelength  of  illuminating  beam 
(6328  A),  and 

V3  =  the  voltage  applied  to  the  electrode. 

From  these  equations  it  can  be  seen  that  the  volt¬ 
age  required  for  2 it  phase  modulation  is  dependent  on  the 
length-to-thickness  ratio  (*2^3)  the  modulator  ele¬ 
ments.  For  the  present  modulator  dimensions,  approxi¬ 
mately  27  volts  is  required  for  a  complete  2 tr  phase  modu¬ 
lation  when  the  light  is  polarized  parallel  to  the  X3  axis, 
and  approximately  120  volts  is  required  when  the  light  is 
polarized  parallel  to  the  Xj  axis.  In  order  to  be  within 
the  voltage  range  of  high-speed  transistor  electronics, 
the  illuminating  polarization  of  the  crystal  modulator 
was  chosen  to  be  parallel  to  the  X3  axis. 
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The  electric  field  distribution  is  nonuniform  and 
also  not  parallel  to  the  X3  axis  at  positions  which  are  not 
near  the  center  area  under  the  electrodes  due  to  the  elec¬ 
tric  field  "fringing"  effects.  Therefore,  light  which  prop¬ 
agates  along  these  parts  of  the  crystal  will  not  be  uni¬ 
formly  phase  modulated  and  also  will  become  elliptically 
polarized.  A  mask  is  used  to  block  the  light  from  all  but 
the  central  areas  under  the  electrodes.  The  present  mask 
consists  of  a  linear  array  of  I00-by-100-^m  apertures 
located  on  the  same  centers  as  the  electrodes  (Fig.  7). 

The  relative  dielectric  constant  (cr)  along  the  X3 
axis  for  lithium  tantalate  is  43  (Ref.  4),  and  therefore 
the  effective  capacitance  of  the  elements  (which  limits 
the  frequency  response  of  the  modulator)  is  approximately 
10  pF. 

Measurements  on  the  masked  modulator  indicate 
that  approximately  32  volts  is  required  for  2rr  phase  modu¬ 
lation  and  that  the  modulation  follows  a  linear  voltage  law. 

With  no  voltages  applied  to  the  crystal  electrodes, 
phase  distortion  of  an  illuminating  plane  wave  has  been 
observed  at  the  output  surface  of  the  crystal.  This  phase 
distortion  (caused  by  nonuniformities  of  the  crystal)  can 
be  compensated  by  applying  appropriate  voltages  to  the 
modulator  electrodes. 

The  present  crystal  modulator  has  46  continguous 
useful  independent  elements.  There  are  several  prob¬ 
lems  associated  with  attempting  to  increase  the  number 
of  elements  in  the  array.  For  2ir  phase  voltages  which 
are  compatible  with  high  speed  transistor  electronics, 
the  length-to -thickness  ratio  of  the  modulator  elements 
(13/^3  in  Eq.  (6))  must  be  large.  The  difficulty  in  fab¬ 
ricating  thin  crystals  has  limited  the  crystal  lengths  to 
approximately  25  mm  (for  0. 1  mm  thickness).  A  higher 
density  of  electrodes  on  the  present  crystal  will  not 
result  in  a  satisfactory  modulator.  If  the  electrodes  are 
located  on  closer  centers,  the  cross  talk  between  elements 
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(the  field  under  one  electrode  due  to  a  voltage  applied  to 
an  adjacent  electrode)  becomes  high.  If  the  width  of  the 
electrode  is  decreased,  the  electric  field  "fringing"  ef¬ 
fects  becomes  a  limiting  factor  because  the  useful  area 
where  the  field  is  uniform  becomes  smaller. 

An  attempt  has  been  made  to  use  two  25-mm 
crystals  mounted  end  to  end  in  order  to  make  a  longer 
array.  However,  the  axes  of  these  crystals  could  not  be 
aligned  with  each  other.  Even  if  the  crystal  axes  could 
be  aligned,  the  physical  length  of  the  array  (50  mm) 
would  require  a  high  quality  transform  lens  in  order  to 
avoid  distortion  in  the  transform  plane, 

A  more  detailed  discussion  of  the  electro -optic 
effect  and  the  derivation  of  Eqs.  (6)  and  (7)  are  given  in 
Appendix  C. 
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BEAM-STEERING  PROPERTIES  OF  AN  OPTICAL 
PHASED-ARRAY  ANTENNA  * 


For  many  years,  phasedrarray  antennas  have  been 
used  in  radar  to  perform  two-dimensional  electronic  beam 
steering  at  microwave  frequencies.  With  th^  development 
of  optical  phased  arrays,  this  electronically  controlled 
beam-steering  technique  can  now  be  extended  to  optical 
wavelengths.  In  this  section,  the  beam-steering  proper¬ 
ties  of  an  optical  phased  array  will  be  briefly  described. 

A  more  detailed  description  is  given  in  Appendix  D. 

I 

Consider  first  a  linear  array  of  coherent  iso¬ 
tropic  light  sources  in  the  front  focal  plane  of  ^  "perfect" 
lens  (Fig.  8).  The  spatial  Fourier  transform  (or  an* 
tenna  pattern)  of  this  array  exists  ip  the  back  focal  plane 
of  the  lens  and  is  described  by  the  equation: 


Ua(£,t?)  =  £  Am  exp(;j<pm)  exp  -j  ~(m  -  -jj— )  AX$J,  (B); 


where  U  (£,?))  =  the  complex  light  distribution  in  the 
d  Fourier  transform  plane  (note  that 

U  is  independent  of  7)  for  this  one-: 
dimensional  array), 

M  =  the  number  of  point  radiators, 

i  th  1 

A  =  the  amplitude  of  the  m  radiator, 

m  1  th 

o  =  the  relative  phase  of  the  m 

m 

radiator, 

AX  =  the  centen  spacing  of  the  elements, 

X  =  the  wavelength  of  the  illumination, 

and 

f  =*  the  focal  length  of  the  lens. 


-  19  - 


TMK  JOHNS  HOMUMI  UNIVUMTV 

APPLIED  PHYSICS  LABORATORY 

•HVta  >MIII>.  MMmjWMI 


PiB.  8  LINEAR  ARRAY  OF  M  EQUALLY  SPACED  ISOTROPIC  RADIATORS  IN  THE 
FRONT  FOCAL  PLANE  OF  A  LENS 


T»*  JOHN*  HOMIM  JN1V*1**JTY 

APPLIED  PHYSICS  LABORATORY 

**m»  *mn.  MAorvut* 


When  the  radiators  are  of  the  same  amplitude 
and  phase,  this  equation  reduces  to: 


ua(«> =  A 


sin  [(MirAX£)/Xj»] 
sin  [(ffAX5)/X/l  ‘ 


(9) 


The  power  in  the  transform  plane  is  given  by: 


I  ua(€>l 2  =ua($)ua*<s> 

_  ^2  sin2  [(MffAX£)/X/]  (1( 

sin2  [(ffAX£)/X/] 

A  plot  of  Eq.  (10)  for  six  point  radiators  (M  =  6) 
is  shown  in  Fig.  9.  The  grating  lobes  (or  ambiguity 
peaks)  of  the  pattern  occur  when  both  the  numerator  and 
denominator  vanish  simultaneously.  This  occurs  when¬ 
ever: 


(ID 


where  K  =  0,±  1,  ±2,  ...  These  peaks  represent  the 
ambiguities  of  the  transform  pattern. 


The  relative  power  at  the  peaks  is  proportional 
to  M^A^,  and  their  width  (from  the  peak  to  the  first  zero) 
is: 


V 

MAX  * 


(12) 


Consider  the  case  where  a  linear  phase  slope  of 
p  cycles  exists  across  the  array.  The  phase  of  the  mth 
element  is  defined  by 


(13) 
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When  this  is  substituted  into  Eq.  (8),  the  result  is 

”.«»  -  *  sg^r  •  (i*: 

This  equation  has  a  peak  whenever: 

«  =  <K+^-  <15 

where  K  =  0,  ±1,  ±2,  ... 

This  pattern  has  the  same  shape  as  Fig.  9  except 

that  now  the  peaks  are  shifted  along  the  £  axis  by  p(^£^). 

Therefore,  the  location  of  the  peaks  on  the  £  axis  is  pro¬ 
portional  to  the  phase  slope  across  the  array.  Using  the 
Rayleigh  resolution  criteria  (that  the  peak  of  one  beam  be 
in  the  null  of  an  adjacent  beam),  there  are  M  resolvable 
beam  positions  between  the  ambiguities. 

A  similar  analysis  can  be  used  to  determine  the 
pattern  of  a  two-dimensional  array  of  isotropic  radiators. 
Consider  an  M-by-N  array  of  isotopic  radiators  where 
the  spacing  of  the  elements  in  the  X  direction  is  AX  and 
the  spacing  of  the  elements  in  the  Y  direction  is  AY 
(Fig.  10).  The  two-dimensional  array  pattern  has  the 
form: 


Ua(C.tf 


E  Z  A  exp(j(Pmn)  exp(-j~«m  +(n AY?j}).(16) 

m*l  n=l  J 


When  the  radiators  are  of  the  same  amplitude  and 
phase,  the  power  in  the  transform  plane  is  described  by 


tM^r>>|2  =  A2 


sin  |(MirAX£)/X/]  2  sin  [(Nff AYt?)/X/]  2 
sin  (UrAX£ j/X/3  j  [sin  [(irAYTfj/X/ ] 


(17) 
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Fig.  10  TWO-DI MENSIONAL  OPTICAL  PHASED  ARRAY 
IN  THE  FRONT  FOCAL  PLANE  OF  A  LENS 


TM«  JOHN*  MOPKM*  UNtVflWTV 

APPLIED  PHYSICS  LABORATORY 
■nvn  IMM,  MUIUM) 


Figure  11  shows  the  width  {from  null  to  null)  and  spacing 
of  the  grating  lobes.  {For  simplicity,  the  sidelobes  have 
not  been  included. )  The  widths  of  the  lobes  (from  the 
peak  to  the  first  null)  in  the  £  and  T)  directions  are: 


Kf 

MAX 


and 


A?) 


X/ 

NAY  * 


(18) 


The  relative  power  of  the  peaks  is  proportional  to 


Consider  now  a  linear  two-dimensional  phase  slope 
across  the  array  where  the  phase  of  the  (m,  n)  element  is 
described  by 


<P 


mn 


2trm  2?rn 

p  +  _rT 


q . 


(19) 


The  plane  of  constant  phase  is  tilted  such  that  there  are 
p  cycles  of  phase  across  the  array  in  the  X  direction 
and  q  cycles  of  phase  across  the  array  in  the  Y  direc¬ 
tion.  Substituting  Eq.  (19)  into  Eq.  (16),  the  resulting 
power  in  the  transform  plane  is  described  by: 


f  sin  [(MitAX£)/X/ -  irp] 

2 

sin  [(NttAYt})/\/  -  irq] 

[sin  [0rAX$)/x/  -  rrp/Mj 

sin  [(trAY -  ffq/N) 

2 

(20) 


The  pattern  takes  the  same  shape  as  Fig.  1 1  except  that 

the  peaks  are  now  shifted  along  the  £  axis  by  ) 

and  along  the  T)  axis  by  q(-^^).  Therefore,  the  location 

of  the  peaks  in  the  (£,?))  plane  indicates  the  two-dimen¬ 
sional  phase  slope  across  the  array. 
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Again  using  the  Rayleigh  resolution  criteria,  there 
are  M  resolvable  beam  positions  in  the  £  direction  and 
N  resolvable  beam  positions  in  the  rj  direction.  This 
results  in  a  total  of  M  times  N  resolution  cells  in  the 
two-dimensional  area  between  the  grating  lobes. 
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4.  ELECTRONICS 


The  electronics  required  for  performing  optical 
beam  steering  with  the  Membrane  Light  Modulator  (and 
later,  with  a  slight  modification,  the  crystal  modulator) 
will  be  described  in  this  section.  The  use  of  these  elec¬ 
tronics  for  performing  radar  signal  processing  will  be 
discussed  in  Section  8. 

In  order  to  perform  beam  steering  with  a  phased - 
array  antenna,  the  appropriate  linear  phase  slope  must 
be  applied  across  the  array  (as  discussed  in  Section  3). 
The  MLM  electronics  shown  in  Fig.  12  generates  the 
voltages  for  each  row  of  the  MLM  that  is  required  to 
steer  the  beam.  For  a  given  beam  position,  the  appro¬ 
priate  digital  voltage  value  for  each  row  of  the  MLM  is 
generated  by  either  the  Beam  Position  Generator  (BPG) 
or  by  a  digital  computer.  These  voltage  values  are  accu¬ 
rate  to  6  bits.  Paper  tape  can  be  used  to  store  the  com¬ 
puter-generated  voltage  values. 

Two  levels  of  registers  (input  and  storage)  are 
used.  One  set  of  data  is  held  in  the  storage  registers 
(one  6-bit  register  per  MLM  row)  while  a  new  set  of  data 
is  loaded  into  the  input  registers.  When  the  input  regis¬ 
ters  are  completely  loaded,  the  new  set  of  data  is  trans¬ 
ferred  in  parallel  into  the  storage  registers.  The  out¬ 
puts  of  the  storage  registers  are  directly  connected  to 
the  MLM  electrodes  via  the  D/A  converters  and  ampli¬ 
fiers. 


The  computer -gene  rated  words  (and  also  the 
paper -tape  words)  are  loaded  sequentially  into  the  input 
registers.  The  BPG-generated  words  can  be  loaded 
sequentially.  To  provide  for  fast  beam  stepping,  the 
BPG  generated  words  can  also  be  loaded  into  the  input 
registers  in  a  quasi -parallel  mode  that  consists  of 
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Fig.  12  BLOCK  DIAGRAM  OF  THE  HYBRID  CORRELATOR  ELECTRONICS 
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parallel  groups  of  8  channels  (with  4  channels  remaining 
for  the  last  group).  The  8  channels  within  a  group  are 
loaded  in  sequence.  In  this  manner,  all  100  channels 
are  loaded  with  8  clock  pulses  rather  than  the  100  clock 
pulses  required  in  a  purely  sequential  operation. 

Since  the  phase  of  the  light  reflected  by  the  MLM 
elements  is  proportional  to  the  square  of  the  voltages 
applied  to  the  elements,  the  BPG  and  the  computer  must 
first  generate  the  phase  slope  required  for  a  particular 
beam  position,  and  then  take  the  square  root  of  these 
values  in  order  to  obtain  the  relative  voltages  required. 
The  phase  slopes  are  generated  modulo  2 ir  phase. 

With  100  elements  in  the  MLM,  there  are  100  re¬ 
solvable  beam  positions  between  ambiguities  (using  the 
Rayleigh  criterion).  The  computer  can  be  used  to  gen¬ 
erate  the  voltages  required  for  these  100  beam  positions. 
The  computer  is  also  used  to  generate  the  voltages  for 
two-dimensional  beam  steering  with  the  thinned  MLM 
arrays  (discussed  in  Section  6).  The  computer  is  also 
used  to  determine  the  phases  and  MLM  voltages  of  the 
radar  pulses  (discussed  in  Section  8).  The  computer  out. 
put  can  either  be  loaded  directly  into  the  MLM  input 
registers  or  stored  on  paper  tape. 

The  number  100  is  difficult  to  work  with  in  a 
digital  simulator,  so  for  the  purpose  of  the  BPG  elec¬ 
tronics,  the  unambiguous  region  between  grating  lobes 
is  divided  into  128  beam  positions.  These  128  positions 
are  defined  as  the  unit  beam  positions  and  can  be  further 
subdivided  into  fractional  positions  of  1/2,  1/4,  and  1/8. 

The  BPG  has  several  modes  of  operation.  In  the 
free -running  mode,  the  MLM  is  stepped  continuously 
through  all  the  beam  positions  in  unit  or  fractional  (1/2, 
1/4,  or  1/8)  steps.  The  stepping  speed  is  controlled  by 
the  basic  clock  frequency.  The  minimum  time  per  step 
is  limited  by  the  digital  logic  to  approximately  2  j*s.  In 
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the  selected  beam  position  mode,  the  simulator  free  runs 
until  it  reaches  the  selected  beam  position  and  then  stops. 
Provisions  have  been  made  to  allow  digital  noise  to  be 
added  from  an  external  source  to  the  BPG- generated  phase 
words. 


The  computer  is  used  to  generate  the  voltage  values 
required  for  beam  steering  with  the  multi-channel  lithium 
tantalate  phase  modulator.  A  slight  modification  of  the 
output  amplifiers  is  required  in  order  to  compensate  for 
the  different  2 tr  voltages  of  the  MLM  (84  volts)  and  the 
crystal  (32  volts).  The  voltages  required  to  compensate 
for  the  phase  distortion  across  the  crystal  (discussed  in 
Section  2)  are  included  in  these  computer  calculations. 
Although  the  BPG  has  not  been  used  for  beam  steering  with 
the  crystal,  it  could  be  used  if  the  square  root  operation 
(which  is  required  for  the  MLM)  were  removed  and  the 
appropriate  phase  compensation  offsets  added. 

A  detailed  description  of  the  simulator  electronics 
is  given  in  Appendix  B. 
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5.  ONE-DIMENSIONAL  OPTICAL  BEAM  STEERING 


Both  the  Membrane  Light  Modulator  ana  the  multi¬ 
channel  crystal  modulator  have  been  used  to  perform  one¬ 
dimensional  phased-array  optical  beam  steering.  The 
optical  configuration  for  these  two  modulators  are  dif¬ 
ferent  and  will  be  described  in  this  section.  The  one¬ 
dimensional  beam  steering  results  for  the  two  modula¬ 
tors  will  also  be  given. 

MLM  BEAM  STEERING 

A  schematic  of  the  optical  configuration  for  one¬ 
dimensional  beam  steering  with  the  MLM  is  shown  in 
Fig.  13.  A  He-Ne  laser  (6328  A)  is  used  as  the  coherent 
light  source.  Conventional  collimating  optics  expand 
the  laser  beam  and  illuminate  the  MLM.  A  beam  splitter 
is  used  to  direct  the  light.  The  light  which  reflects  off 
the  MLM  is  imaged  via  the  beam  splitter  and  lens  onto 
the  transmit  mask  at  unity  magnification.  (The  purpose 
of  the  transmit  mask  was  discussed  in  Section  2. )  The 
alignment  of  the  transmit  mask  and  the  modulator  image 
is  very  critical.  A  slight  error  in  magnification  of  the 
image  or  position  of  the  mask  will  result  in  a  poor  mask¬ 
ing  of  the  modulator  elements.  As  described  in  Appendix 
E,  unity  magnification  is  the  best  choice  for  simplifica¬ 
tion  of  the  alignment  procedure. 

The  image  of  the  masked  MLM  exists  in  the  front 
focal  plane  of  the  transform  lens,  and  the  spatial  Fourier 
transform  of  the  light  transmitted  through  the  mask  ex¬ 
ists  in  the  back  focal  plane  (output  plane). 

The  light  intensity  in  the  output  plane  is  read  by 
a  scanning  mirror  readout  system.  The  mirror  scans 
a  projected  image  of  the  transform  pattern  across  a  slit 
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in  front  of  the  photomultiplier  so  that  the  current  output 
of  the  photomultiplier  as  a  function  of  time  corresponds 
to  the  intensity  of  the  light  in  the  transform  pattern  as  a 
function  of  the  position  along  the  scan.  The  photomulti¬ 
plier  output  is  displayed  on  an  oscilloscope. 

Since  all  the  elements  along  one  row  of  the  MLM 
are  controlled  by  a  single  electrode,  linear  phase  slopes 
(using  the  full  array)  can  be  obtained  only  in  a  direction 
perpendicular  to  the  electrodes  (which  corresponds  to  the 
X  direction  in  the  input  plane  since  the  electrodes  have 
been  aligned  to  be  parallel  to  the  Y  axis).  Therefore 
one-dimensional  beam  steering  occurs  along  the  £  axis 
ip  the  transform  plane.  The  masked  MLM  elements  have 
a  finite  size,  and  therefore  Eq.  (20)  (which  describes  the 
transform  pattern  for  an  array  of  point  radiators)  must 
be  multiplied  by  the  element  pattern.  The  transform  of 
a  uniformly  illuminated  circular  aperture  is  the  familiar 
Airy  pattern,  and  therefore  the  power  in  the  transform 
plane  for  the  masked  MLM  array  can  be  approximated  by: 


i 


sin  ((MffAX£)/X/  -  up] 
sin  ((ffAX£)/X/  -  ffp/MJ 


2 

sin  [M ffAXp/X/1 

sin  [ffAXrfTxT)  '^U 


where  R  =  the  radius  of  the  mask  apertures  (20  pm), 

p  =  T7  (£2  + 1)2)^#  and  (for  the  MLM) 

M  =  N  =  100  and 
AX  =  AY  =  50  pm. 

A  photograph  of  the  transform  pattern  for  the 
masked  MLM  is  shown  in  Fig.  14a.  (The  transform  pat¬ 
tern  is  rotated  by  45°  corresponding  to  the  actual  orien¬ 
tation  of  the  MLM  array,  which  is  at  45°.  )  Figure  14b  is 
a  plot  of  the  intensity  along  the  £  axis  (marked  A -A  in 
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Fig.  14  MLM  TRANSFORM  PATTERN 
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the  figure).  From  Eq.  (21),  there  are  seven  grating  lobes 
under  the  main  lobe  of  the  element  pattern.  The  distance 
between  the  grating  lobes  represents  the  distance  the  beam 
can  be  steered  unambiguously. 

A  photograph  of  the  main  grating  lobe  on  an  ex¬ 
panded  scale  is  shown  in  Fig.  15.  A  diode  is  used  as  the 
load  for  the  photomultiplier  in  order  to  approximate  a 
logarithmic  amplifier.  The  theoretical  intensities  of  the 
sidelobes  relative  to  the  main  lobe  can  be  obtained  from 
Eq.  (21).  Table  1  compares  the  theoretical  intensities 
of  the  sidelobes  to  those  measured.  As  can  be  seen,  the 
measured  values  compare  favorably  with  the  theoretical 
values. 


Table  1 

Comparison  of  Theoretical  and  Measured  Sidelobe 
Intensities  for  MLM 


Lobe 

Center  Lobe  (reference) 
First  Sidelobe 
Second  Sidelobe 
Third  Sidelobe 


Theoretical 

Intensity 

(dB) 

0 

-13.2 

-17.7 

-20.8 


Measured 

Intensity 

(dB) 

0 

-13.2 

-17.6 

-20.2 


Due  to  the  frequency  response  of  the  diode -photomulti¬ 
plier  combination  and  also  the  finite  size  of  the  slit  in 
front  of  the  photomultiplier,  the  nulls  in  Fig.  15  are  not 
a  good  representation  of  the  actual  nulls  achieved  in  the 
transform  pattern. 

Since  the  MLM  has  100  independent  rows  (M  =  100), 
there  are  100  Rayleigh  beam  positions  between  the  grating 
lobes.  Figure  16  shows  several  adjacent  Rayleigh  beam 
positions.  Notice  that  the  peak  of  one  beam  occurs  in  the 
first  null  of  the  adjacent  beam  position  (satisfying  the 
Rayleigh  resolution  criteria).  The  intensity  at  the  crossing 
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Fig.  1 5  CENTRAL  GRATING  LOBE  OF  THE  MLM  TRANSFORM  PATTERN 
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Fig.  16  MLM  ONE-DIMENSIONAL  BEAM -STEERING  RESULTS  SHOWING  SEVERAL 
ADJACENT  RAYLEIGH  BEAM  POSITIONS  SUPERIMPOSEO  ON  ONE  PHOTO¬ 
GRAPH 
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point  is  observed  to  be  4  dB  below  the  intensity  of  the 
peak  of  one  of  the  beams  (obeying  Eq.  (21)),  If  a  3-dB 
beamwidth  were  used  as  the  resolution  criteria,  113 
beam  positions  could  be  resolved  between  the  ambiguities. 

A  photograph  of  all  100  beam  positions  is  shown  in 
Fig.  17.  The  averaged  envelope  follows  the  Airy  disk 
element  pattern.  However,  a  large  variation  of  the  peak 
intensity  from  one  beam  position  to  the  next  is  observed. 
(This  is  especially  evident  for  beam  position  number  0. ) 

This  is  because  the  element  pattern  varies  with  voltage. 
When  no  voltage  is  applied  to  the  modulator  elements, 
the  wavefront  illuminating  the  circular  mask  apertures 
is  flat,  and  the  transform  pattern  is  the  Airy  pattern. 
However,  when  a  voltage  is  applied,  the  modulator  ele¬ 
ments  deflect  as  paraboloids.  The  wavefront  illuminat¬ 
ing  the  individual  mask  apertures  now  appears  parabolic. 
This  results  in  a  widening  of  the  main  lobe  of  the  element 
pattern  and  also  a  lowering  of  the  peak  intensity  (along 
boresight).  The  larger  the  applied  voltage,  the  greater 
the  effect  since  the  wavefront  is  more  parabolic.  The 
theoretical  intensity  on  boresight  as  a  function  of  peak 
phase  (for  20  ym  diameter  mask  apertures)  is  plotted 
in  Fig.  18.  By  applying  the  same  voltage  to  all  the  MLM 
electrodes,  the  variation  of  intensity  on  boresight  as  a 
function  of  the  peak  phase  can  be  measured.  These  mea¬ 
surements  are  also  plotted  in  Fig.  18.  Up  to  around 
3ff/4  phase,  the  theoretical  and  empirical  curves  are 
approximately  the  same.  However,  beyond  3jt/4  phase 
the  measured  intensity  on  boresight  drops  very  quickly. 

This  is  attributed  to  abberations  in  the  imaging  lens. 

The  theoretical  plot  in  Fig.  18  was  made  assuming  a 
perfect  lens  of  infinite  aperture  (a  point  impulse  re¬ 
sponse  function  for  the  lens).  For  a  perfect  lens  with  a 
finite  aperture,  the  impulse  response  function  is  a  sym¬ 
metric  Airy  pattern.  For  an  f! 4  lens,  the  theoretical 
intensity  versus  phase  curve  has  been  found  to  follow 
fairly  close  to  the  infinite  aperture  curve  for  all  phase 
modulations.  The  impulse  response  function  for  the  imaging 
lens  in  the  system  (which  is  // 4.  8)  has  been  measured  and 
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NOTE;  FCR  PHOTOGRAPHIC  REASONS,  THE  0  LEVEL  HAS  BEEN 
SHIFTED  1  CM  BELOW  THE  BOTTOM  SCALE  LINE.  THE 
AMBIGUITY  SPACING  CORRESPONDS  TO  B  CM.  THE 
BORESIGHT  BEAM  (AT  £  =  0)  HAS  BEEN  CENTERED  ON 
THE  PHOTOGRAPH 


Fig.  17  THE  100  RAYLEIGH  BEAM  POSITIONS  OF  THE  MLM 
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Fig.  18  INTENSITY  AT  BORESIGHT  OF  THE  MLM  TRANSFORM  PATTERN  AS  A 
FUNCTION  OF  THE  AVERAGE  PHASE  ACROSS  THE  MLM  ELEMENTS 
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is  shown  in  Fig.  19.  The  fact  that  this  is  an  imperfect  lens 
(as  evident  from  the  asymmetrical  impulse  response  func¬ 
tion)  results  in  a  degradation  of  the  image  and  probably  is 
responsible  for  the  discrepancy  in  the  measured  intensity 
modulation  above  3rr / 4  phase. 


CRYSTAL  BEAM  STEERING 

A  schematic  of  the  optical  configuration  for  one¬ 
dimensional  beam  steering  with  the  multi-channel  crystal 
modulator  is  shown  in  Fig.  20.  A  He-Ne  laser  (6328  A) 
is  used  as  the  coherent  light  source.  A  polarizer  is  used 
to  cancel  out  the  depolarizing  effects  of  the  quartz-coated 
mirrors.  The  polarization  rotator  is  used  to  align  the 
polarization  of  the  input  beam  parallel  to  the  X3  axis  of 
the  crystal. 

Conventional  spherical  collimating  optics  expand 
the  beam  to  illuminate  the  crystal.  Because  of  the  rec¬ 
tangular  cross  section  of  the  crystal  (0.  1  by  25  mm),  only 
a  small  percentage  of  the  circular  beam  is  actually  used. 
By  using  cylindrical  optics,  the  circular  beam  can  be 
compressed  in  one  dimension  such  that  the  beam  cross 
section  more  closely  corresponds  to  the  crystal  cross 
section  (Fig.  21).  The  amount  of  compression  is  pro¬ 
portional  to  the  ratio  of  the  effective  focal  lengths  of  the 
cylindrical  lenses: 


m 


{_2 
*1  ’ 


(22) 


As  discussed  in  Section  2,  a  mask  is  required  to 
block  the  light  from  the  areas  between  the  modulator  elec¬ 
trodes.  The  present  mask  consists  of  a  linear  array  of 
lOO-by-lOO-ym  apertures  which  are  on  500-ym  centers 
(Fig.  22).  The  transmit  mask  is  positioned  within  100  ym 
of  the  output  surface  of  the  crystal  by  means  of  a  3-axis 
translator. 
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Fig.  19  MEASURED  IMPULSE  RESPONSl 
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Fig.  22  TRANSMIT  MASK  FOR  46  CHANNEL  CRYSTAL  PHASE  MODULATOR 
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The  output  surface  of  the  masked  crystal  is  located 
in  the  front  focal  plane  of  the  transform  lens,  and  the 
spatial  Fourier  transform  of  the  array  exists  in  the  back 
focal  plane  (transform  plane).  The  light  intensity  in  the 
transform  pattern  is  read  by  a  scanning  mirror  readout 
system.  The  mirror  scans  a  projected  image  of  the  trans¬ 
form  past  a  slit  in  front  of  the  photomultiplier  so  that  the 
current  out  of  the  photomultiplier  as  a  function  of  time 
corresponds  to- the  light  intensity  in  the  transform  plane 
as  a  function  of  the  distance  along  the  scan.  The  photo¬ 
multiplier  output  is  displayed  on  an  oscilloscope. 

The  multi-channel  crystal  modulator  can  be  con¬ 
sidered  as  a  linear  phased -array  antenna.  The  transform 
pattern  of  a  linear  array  of  isotropic  radiators  was  given 
in  Eq.  (14).  Since  the  crystal  elements  are  rectangular, 
Eq.  (14)  must  be  multiplied  by  the  element  transform  pat¬ 
tern.  The  power  in  the  transform  plane  of  the  multi- , 
channel  crystal  modulator  is  described  by  the  equation: 


where  a  =  the  width  of  the  crystal  elements  in  the  X 
direction,  and 

b  =  the  width  of' the  crystal  elements  in  the  Y 

direction  (for  the  .present  mask,  a  =  b  =  100  4m). 

*  1 

For  the  present  crystal,  there  are  46  useful  ele¬ 
ments  (M  =  46),  and  therefore  there  are  46  Ilayleigh  beam 
positions  between  the  grating  lobes.  A  photograph  of  sev¬ 
eral  adjacent  beam  positions  is  shown  in  Fig.  23.  v  The 
width  of  the  lobes  obey  Eq.  (12)  and  the  ihtensity  at  the 
crossing  point  between  two  adjacent  beam  positions 
(-4  dB)  satisfies  Eq.  (23).  :  ? 


48 


1 


^  Vg  j1^  WBWfMgW>.lffCWt  J *  %•;  kWWiJj.Iy'riWM*^WW1  UJ!.L.I^WWW^WWWW 


’AWMMmw 


«nnmm»')f*  v*.j*ii'ijmm\mKm'»K*c-**iKa 


Tmc  jphnb  Menuwi  uwivmtirr 
APPLIED  PHYSICS  LABORATORY 

Sm.v*»  >n»iw  kUimjM<« 


$/(V/48AX) 

Fig.  23  ONE-DIMENSIONAL  BEAM-STEE RING  RESULTS 

USING  THE  MULTICHANNEL  CRYSTAL  MODULATOR 
{SEVERAL  ADJACENT  RAYLEIGH  BEAM  POSITIONS 
ARE  SUPERIMPOSED  ON  ONE  PHOTOGRAPH) 


A  photograph  of  all  46  beam  positions  is  shown  in 
Fig.  24.  Unlike  the  MLM  (Fig.  17),  the  intensity  is  fairly 
uniform  for  all  of  the  beam  positions.  This  is  because 
the 'amplitude  of  the  crystal  elements  is  not  a  function  of 
the  applied  phase  modulation.  In  other  words,  the  crystal 
acts  as  a  piston  modulator. 

A  photograph  of  the  grating  lobe  structure  is  shown 
in  Fig.  25.  From  Eq.  (23),  there  are  nine  grating  lobes 
under  the  main  lobe  of  the  element  pattern. 
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Fig.  24  RAYLEIGH  BEAM  POSITIONS  OVER  ENTIRE  AMBIGUITY 

INTERVAL  FOR  THE  46  CHANNEL  CRYSTAL  PHASE  MODULATOR 
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Fig.  25  GRATING  LOBE  STRUCTURE  OF  THE  CRYSTAL  MODULATOR 
TRANSFORM  PATTERN 
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6.  TWO-DIMENSIONAL  OPTICAL  BEAM  STEERING 


As  discussed  in  Section  3,  the  requirement  for  two- 
dimensional  beam  steering  with  a  phased-array  antenna  is 
the  capability  of  generating  two-dimensional  phase  slopes 
across  the  antenna.  When  an  appropriate  optical  system 
is  employed,  both  the  MLM  and  crystal  modulators  can 
be  used  to  perform  two-dimensional  beam  steering  (in  this 
case,  only  N  elements  are  used  from  the  N-by-N  array). 
Two-dimensional  beam  steering  has  been  successfully  per¬ 
formed  using  the  MLM.  A  description  of  the  MLM  system 
and  results  will  be  given  in  this  section.  Also,  an  optical 
system  for  performing  two-dimensional  beam  steering  with 
the  crystal  modulator  will  be  described. 


OPTICAL  CONFIGURATION  FOR  TWO-DIMENSIONAL 
BEAM  STEERING  WITH  THE  MLM 

Since  the  elements  along  each  row  of  the  MLM 
array  are  controlled  by  a  single  electrode,  the  phase  of 
each  of  these  elements  is  the  same.  Phase  slopes  parallel 
to  the  electrode  structure  are  required  in  order  to  perform 
two-dimensional  beam  steering,  and  therefore  only  one 
element  from  each  row  can  be  used.  A  mask  is  used  to 
block  the  light  from  all  but  the  appropriate  element.  A 
typical  mask  is  shown  in  Fig.  26a.  For  simplicity,  a 
20-by-20  array  has  been  shown.  For  this  example,  the 
light  is  blocked  from  all  but  the  first  element  in  the  first 
row,  the  eleventh  element  in  the  second  row,  etc.  This 
thinned  array  can  be  used  to  obtain  a  two-dimensional 
tilt  of  the  plane  of  equal  phase  as  shown  in  Fig.  26b.  The 
optical  configuration  for  two  dimensional  beam  steering 
with  the  MLM  (Fig.  26c)  is  the  same  as  for  one-dimen¬ 
sional  beam  steering  (Fig.  13)  except  that  now  the  thinned 
array  transmit  mask  is  used.  When  the  appropriate 
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(a)  TRANSMIT  MASK 


(b>  TWO-DIMENSIONAL  PHASE  SLOPE 
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Fig.  26  TWO  DIMENSIONAL  OPTICAL  BEAM  STEERING  USING  THE  MLM 


-  54  - 


THC  JOHNS  HOSKINS  UNIVUSITY 

APPLIED  PHYSICS  LABORATORY 

SKVtM  SNMt.  tUITUlM 


voltages  are  applied  to  the  modulator  electrodes,  the  beam 
from  the  thinned  MLM  array  can  be  steered  in  two  dimen¬ 
sions. 


SELECTION  OF  THE  ELEMENT  LOCATIONS  IN 
THE  THINNED  MLM  ARRAYS 

If  independent  access  could  be  made  to  each  ele¬ 
ment  in  the  100-by-l00  MLM  array,  the  beam  could  be 
steered  into  10  000  Rayleigh  resolvable  beam  positions. 
When  the  locations  of  the  thinned  array  elements  are 
chosen  judiciously,  this  resolution  of  the  whole  array  is 
maintained.  One  of  the  criteria  used  for  selecting  the 
best  location  of  the  elements  in  the  chinned  MLM  array 
was  that  the  resolution  in  the  transform  pattern  of  the 
entire  100-by-lU0  array  should  be  maintained. 

In  order  to  avoid  the  observation  of  false  beam 
positions,  it  is  necessary  to  minimize  the  peak  power 
in  lobes  other  than  the  main  grating  lobes.  The  second 
criterion  used  for  selecting  the  best  location  of  the  ele¬ 
ments  in  the  thinned  array  was  that  the  sidelob?  levels 
in  the  transform  pattern  should  be  as  low  as  possible. 

For  a  perfect  "thumbtack”  response  (i.  e. ,  the  sidelobe 
power  is  uniformly  spread  over  the  two-dimensional 
ambiguity  area),  it  can  be  shown  that  the  average  power 
in  the  sidelobes  is  approximately  1/N  times  the  peak 
power  in  the  main  beam.  For  the  thinned  MLM  array, 

N  =  100  and  therefore  the  average  sidelobe  level  is 
20  dB  down  from  the  main  lobe  peak.  However,  it  is 
difficult  to  achieve  a  thumbtack  response  due  to  the 
periodicity  of  the  array  locations  and  the  relatively  small 
number  of  elements. 

The  lower  limit  of  the  peak  sidelobe  level  for  the 
thinned  MLM  arrays  is  determined  by  the  requirement 
that  only  one  element  occur  on  each  electrode.  There¬ 
fore,  the  transform  pattern  along  the  £  axis  is  de- 
scribed  by  a  (sin  £)/£  distribution  in  the  vicinity  of  the 
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main  grating  lobe.  The  peak  intensity  of  the  first  side- 
lobe  for  this  distribution  is  13  dB  down  from  the  peak  in¬ 
tensity  of  the  main  lobe.  The  goal  therefore  is  to  keep 
the  peak  intensity  at  any  other  point  in  the  area  between 
the  grating  lobes  below  this  13-dB  level. 

Equation  (16)  (which  describes  the  transform  pat¬ 
tern  for  a  full  two-dimensional  array  of  point  radiators) 
must  be  modified  to  the  following  form  in  order  to  de¬ 
scribe  a  thinned  array  of  isotropic  radiators  (ignoring 
a  constant  phase  term): 


M_1  /  \ 

U  (4,7})  =  £  A  exp(jtf  )  exp(-jk(mAX4  +  n(m)AY?}]),  (24) 

a  m=0  m  m  '  ' 


where  n(m) 


=  the  relative  column  location  of  the 
radiator  in  the  mth  row, 

=  the  amplitude  of  the  radiator  in  the  mth 
row, 

-  the  phase  of  the  radiator  in  the  mth  row, 
and 


k  =  2 


The  transform  pattern  for  the  boresight  beam  of  this 
thinned  array  can  be  found  by  letting  the  phase  and  am¬ 
plitude  of  each  radiator  be  equal  (also  AX  =  AY  for  the 
MLM): 


M-l 

u  (4,r>)  =  £  exp  l-jk  AX  (m4  +  n(m)r})4  (25) 

a  m=0 


Equation  (25)  can  be  solved  in  closed  form  for  only  a  few 
simple  arrangements  of  the  input  radiators.  Using  Eq. 
(25),  a  computer  program  has  been  written  to  calculate 
the  power  in  the  transform  pattern  for  an  arbitrary 
arrangement  of  the  radiators  in  the  thinned  array.  From 
the  results  of  this  program,  the  "best"  orientation  of  the 
radiators  can  be  chosen. 
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With  a  rather  limited  effort,  no  selection  technique 
has  been  found  that  is  better  than  pseudo -random  orienta¬ 
tions  of  the  elements  in  the  thinned  array.  The  resolution 
of  the  entire  array  is  maintained  and  the  peak  sidelobes 
are  about  10  dB  down  from  the  peak  intensity  of  the  main 
lobe.  Optical  transmit  masks  for  three  different  pseudo¬ 
random  orientations  have  been  obtained. 

A  description  of  the  technique  used  for  selecting 
the  element  locations  for  the  thinned  MLM  arrays  is  given 
in  Appendix  F. 


TWO-DIMENSIONAL  BEAM-STEERING  RESULTS 
USING  THE  MLM 

A  computer  program  has  been  written  to  generate 
the  voltages  required  for  steering  the  beam  of  the  thinned 
MLM  array  to  any  place  in  the  unambiguous  transform 
area.  The  row  and  column  locations  of  the  thinned  array 
elements  and  the  desired  locations  of  the  MLM  beam  are 
read  into  the  computer  program.  The  output  from  the 
computer  program  is  punched  on  paper  tape  in  a  format 
that  is  compatible  with  the  paper  tape  reader  of  the  Digital- 
Optical  Hybrid  Correlator. 

Two-dimensional  optical  beam  steering  has  been 
performed  using  the  three  pseudo -random  optical  masks. 
Photographs  of  the  transform  patterns  for  the  three  masks 
are  shown  in  Fig.  27.  The  transform  pattern  is  at  the 
usual  45°,  corresponding  to  the  MLM  array  orientation 
of  45°.  The  main  grating  lobe  (the  grating  lobe  at  the  cen¬ 
ter  of  the  element  pattern)  has  been  centered  on  these 
photographs. 

Due  to  the  low  light  levels,  a  measurement  of  the 
actual  sidelobe  levels  in  the  transform  plane  could  not 
be  made.  However,  several  observations  can  be  made 
from  the  photographs  of  the  transform  pattern.  First, 
the  light  intensity  in  the  area  between  the  grating  lobes  is 
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(c)  R1  MASK 

Fig.  27  TRANSFORM  PATTERN  FOR  THE  THREE  PSEUDO-RANDOM 
TRANSMITMASKS.  (Cont'd) 
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fairly  uniform  and  therefore  approaches  a  thumbtack  re¬ 
sponse.  Also,  the  sidelobe  patterns  in  the  computer  pro¬ 
gram  plots  (see  Appendix  F)  compare  very  closely  with 
the  sidelobe  patterns  in  the  photographs.  A  dark  line 
structure  connecting  the  grating  lobes  is  also  evident  in 
these  photographs.  This  line  structure  corresponds  to 
the  £  and  t]  axes  of  the  transform  pattern  (and  its  am¬ 
biguities)  along  which  the  transform  pattern  is  described 
(in  the  vicinity  of  the  grating  lobes)  by  a  (sin  x)/x  function. 
The  first  sidelobe  of  this  distribution  can  be  seen  adjacent 
to  the  main  grating  lobes.  Since  the  higher  order  side- 
lobes  of  this  distribution  are  very  low  with  respect  to  the 
light  level  in  the  area  between  the  grating  lobes,  a  dark 
line  structure  will  be  observed  connecting  the  grating  lobes 
in  the  £  and  17  directions. 

Typical  beam  steering  results  for  one  of  the  masks 
(Mask  RRC1)  are  shown  in  Fig.  28.  For  this  example, 
the  beam  was  steered  along  the  diagonal  between  the  grat¬ 
ing  lobes.  A  number  of  comparisons  between  the  steered 
beam  and  the  boresight  beam  (when  zero  phase  slope  ex¬ 
ists  across  the  modulator)  can  be  made: 

1.  The  overall  intensity  in  the  unambiguous 
transform  area  of  the  steered  beam  is  less 
than  that  of  the  boresight  beam. 

2.  The  intensity  of  the  main  steered  lobe  is 
lower  than  the  intensity  of  the  main  boresight 
lobe. 

3.  The  peak  signal  to  peak  sidelobe  power  ratio 
is  lower  in  the  steered  beam. 

4.  The  transform  pattern  of  the  steered  beam 
does  not  look  the  same  as  the  transform  pat¬ 
tern  of  the  boresight  beam.  Sidelobes  occur 
at  different  places  and  the  dark  lines  con¬ 
necting  the  ambiguities  are  no  longer  as  dis¬ 
tinct. 
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(a)  BORESIGHTBEAM  POSITION 
“  <0, 0) 


STEERED  BEAM  POSITION 

N0TE:  MNTHfS 

an^or  the  su"ent  photographs  the  camera  position  was  helo  m«D. 


Fig.  28 


TWO-DIMENSIONAL  BEAM-STEERING  RESULTS  USING  a  THINNED  MLM  ARRAY 
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(c)  STEERED  BEAM  POSITION 


(d)  STEERED  BEAM  POSITION 
lfc  -  (25,25) 


Fig.  28  TWO-DIMENSIONAL  BEAM  STEERING  RESULTS 
USING  A  THINNED  MLM  ARRAY  (cont'd) 
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(•)  STEERED  BEAM  POSITION 
(£rji-<49;49) 


{ft  STEERED  BEAM  POSITION 
(|.  T?) -(99,99) 


Fig.  28  TWO-DIMENSIONAL  BEAM  STEERING  RESULTS 
USING  A  THINNED  MLM  ARRAY  (confd) 
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5.  The  main  lobe  of  the  steered  beam  does  not 
appear  as  distinct  as  the  boresight  lobe. 

; 

6.  A  residual  boresight  lobe  is  not  noticeable 
when  the  beam  is  steered. 

The  observations  listed  as  1  through  5  are  a  result 
of  the  decreasing  amplitude  in  the  element  pattern  as 
higher  voltages  (phases)  are  applied  to  the  MLM  elements 
(discussed  in  Section  5).  For  a  piston  modulator  (without 
this. amplitude  versus  phase  dependence),  the  transform 
pattern  of  the  steered  beam  would  be  essentially  identical 
to  the  boresight  transform  pattern.  The  observation  listed 
as  6  above  is  an  indication  of  a  good  alignment  of  the  trans¬ 
mit  mask. 


TWO-DIMENSIONAL  BEAM  STEERING  USING  THE 
CRYSTAL  MODULATOR 

By  using  appropriate  optics,  a  linear  array  of 
square  radiators  (corresponding  to  the  masked  crystal 
modulator)  can  be  converted  into  a  two-dimensional  array 
as  illustrated  in  Fig.  29.  The  array  in  Fig.  29a  is  de- 
magnified  in  the  X  direction  and  magnitied  in  the  Y 
direction  to  produce  a  compressed  array  of  rectangular 
elements.  As  shown  in  Fig.  29b,  this  new  array  can  be 
masked  with  a  two-dimensional  array  of  circular  aper¬ 
tures  in  a  configuration  similar  to  the  MLM.  Now, 
thinned  array  masks  similar  to  the  masks  designed  for 
two-dimensional  beam  steering  with  the  MLM  can  be 
used  to  perform  two-dimensional  beam  steering  with  the 
crystal. 

The  optics  required  to  obtain  this  compressed 
rectangular  array  are  shown  in  Fig.  30.  Cylindrical 
lenses  are  used  to  demagnify  the  array  in  the  X  direc¬ 
tion  by  an  amount: 


‘X  I. 
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<»>  LINEAR  ARRAY  OF  SQUARE  RADIATORS  (b)  COMPRESSED  ARRAY  OF  RECTANGULAR 

RADIATORS  MASKED  BY  TWO-DIMENSIONAL 
ARRAY  OF  CIRCULAR  APERTURES 


Fig.  29  CONVERSION  OF  A  LINEAR  ARRAY  TO  A  TWO-DIMENSIONAL  ARRAY 
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Fig.  30  ADDITIONAL  OPTICS  REQUIRED  TO  PERFORM  TWO-OIMENSIONAL 
BEAM  STEERING  WITH  THE  CRYSTAL  MODULATOR 
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The  natural  divergence  of  the  modulator  elements  results 
in  a  magnification  in  the  Y  direction  of 


d2  2(1.  22  *5-) 

D  **  D 


as 


5X(/l  +/2) 


D‘ 


(27) 


By  choosing  a  proper  set  of  focal  lengths  for  the  cylindri¬ 
cal  lenses  (/,  and  /  2)»  the  resulting  image  is  similar  to 
Fig.  29b. 


Since  the  crystal  elements  act  like  piston  phase 
modulators,  the  amplitude /phase  dependence  that  is  char¬ 
acteristic  of  the  MLM  is  not  present.  Therefore,  the 
transform  pattern  of  the  steered  beam  should  be  essen¬ 
tially  identical  to  the  transform  pattern  of  the  bore  sight 
beam. 
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7.  CORRELATION  PROPERTIES  OF  AN  OPTICAL 
PHASED  ARRAY 


The  correlation  properties  of  an  optical  phased 
array  will  be  discussed  in  terms  of  the  operations  re¬ 
quired  to  correlate  a  pulse  doppler  waveform.  For  a 
pulse  doppler  transmission,  a  set  of  single  frequency 
pulses  is  transmitted  and  the  echo  phases  for  one  range 
bin  (pulse  width)  are  determined.  For  a  constant  velocity 
target  in  a  noise -free  environment,  there  exists  a  linear 
phase  slope  across  the  return  echoes  (corresponding  to 
the  constant  doppler  frequency  shift).  The  function  of 
the  correlator  is  to  determine  the  phase  slope  (and  there¬ 
fore  the  velocity  of  a  target)  when  noise  is  present.  The 
phase  slope  corresponding  to  a  given  velocity  is  assumed. 
The  correlation  of  the  return  echoes  with  this  linear  phase 
slope  is  given  by: 


C(V)  =  £  Am  exp  (j<Pm)  exp  f-j  4ffVmAt  \  .  (27) 
m=l  '  t  ' 

the  amplitude  of  the  m**1  echo, 

the  phase  of  the  mth  echo, 

the  assumed  velocity  of  the  target, 

the  time  between  successive  pulses  (l/PRF), 
and 

Xt  =  the  wavelength  of  the  transmitted  waveform. 

An  electronic  system  performs  this  correlation  in 
a  brute  force  manner  by  first  rotating  the  vectors  of  the 
returned  pulses  by  the  phase  of  the  assumed  linear  slope 
and  then  performing  a  vector  addition  of  all  the  pulses. 
This  is  then  repeated  for  every  velocity  bin  that  is  in¬ 
vestigated. 


where  A 
<P 


m 


m 
V 

At 


Preceding  page  blank 
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Notice  that  Eq.  (27)  is  similar  in  form  to  Eq.  (8) 
which  describes  the  transform  pattern  for  a  linear  phased 
array  antenna  (ignoring  a  uniform  phase  term).  If  AX  is 
associated  w;th  At,  and  £  is  associated  with  V,  then 
Eq.  (8)  can  be  interpreted  as  a  correlation  function.  The 
phase  of  the  mth  echo  is  impressed  on  the  mth  radiator, 
and  the  spatial  spacing  of  the  radiators  corresponds  to  the 
temporal  spacing  of  the  radar  pulses.  The  position  along 
the  4  axis  in  the  transform  plane  corresponds  to  the  ve¬ 
locity  of  the  target.  The  optical  system  not  only  performs 
the  correlation  instantaneously,  but  also  the  correlation  is 
performed  simultaneously  for  all  velocities. 

The  resolution  and  ambiguities  of  the  pulse  doppler 
waveform  are  the  same  as  those  already  calculated  in  Sec¬ 
tion  3  for  the  phased  array. 

This  correlation  property  can  also  be  explained  in 
terms  of  optical  beam  steering.  When  a  linear  phase 
slope  corresponding  to  a  constant  velocity  target  exists 
across  the  phased  array,  the  beam  is  steered  an  amount 
proportional  to  the  target  velocity. 

Ideally,  the  amplitude  (Am)  and  phase  (<Pm)  of  each 
radar  pulse  should  be  included  in  the  correlation  calcula¬ 
tions  (as  shown  in  Eq.  (27)).  However,  as  discussed  in 
Section  2,  the  MLM  and  crystal  modulators  are  pure 
phase  modulators.  Therefore,  the  amplitude  of  the  input 
data  is  hard  limited.  This  results  in  the  appearance  of 
spurious  signals  and  also  a  degraded  signal-to-noise  ratio 
(Ref.  6’.  Several  techniques  for  maintaining  amplitude 
and  phase  information  with  the  MLM  abd  crystal  modulators 
have  been  proposed  but  not  yet  implemented. 
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8.  ONE-DIMENSIONAL  RADAR  SIGNAL 
PROCESSING 


A  digital-optical  hybrid  correlator  which  has  the 
capability  of  processing  pulse  doppler  and  wideband  radar 
signals  has  been  developed  based  on  the  correlation  prop¬ 
erties  of  the  MLM  phased  array.  This  hybrid  system  re¬ 
tains  the  flexibility  of  the  digital  technology  for  signal 
storage  and  registration  and  the  speed  of  the  coherent 
optical  technology  for  performing  the  correlation  arith¬ 
metic.  Use  of  digital  shift  registers  for  signal  borage 
results  in  a  system  capable  of  operation  at  data  rates 
from  DC  through  the  modulator  response  time. 


PULSE  DOPPLER  WAVEFORMS 

The  use  of  the  digital-optical  hybrid  correlator  to 
process  pulse  doppler  waveforms  is  shown  in  Fig.  31.  As 
discussed  in  Section  7,  a  linear  phase  slope  exists  across 
the  return  echoes  for  a  constant  velocity  target  in  a  noise- 
free  environment.  When  the  phases  for  a  set  of  100  pulses 
are  time  ordered  across  the  100  MLM  rows,  the  output 
beam  is  steered  an  amount  that  is  a  function  of  the  target 
velocity.  For  a  real-life  environment,  the  output  pattern 
corresponds  to  the  correlation  of  the  input  waveform  with 
all  the  possible  phase  slopes  (velocities). 

The  optical  configuration  for  pulse  doppler  signal 
processing  is  the  same  as  the  optical  configuration  for 
one-dimensional  beam  steering  that  was  described  in 
Section  5  (Fig.  13).  A  block  diagram  of  the  receiver  for 
window  mode  processing  is  shown  in  Fig.  32.  For  a 
pulse  doppler  transmission,  the  frequency  switch  is  kept 
at  one  frequency  and  a  set  of  100  pulses  is  transmitted. 

The  phases  of  the  echo  pulses  are  calculated  for  one  range 
bin  (pulse  width),  and  then  the  square  roots  of  these  values 
are  determined  (since  the  MLM  phase  obeys  a  square 
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voltage  law).  These  values  are  then  loaded  into  the  input 
registers  and  ordered  according  to  the  received  time. 

After  the  phase  values  for  all  100  pulses  have  been  loaded 
into  the  input  registers,  they  are  transferred  in  parallel 
to  the  output  registers,  which  are  tied  directly  to  the  MLM 
rows  via  D/A  converters  and  amplifiers. 

The  hybrid  correlator  has  been  used  off-line  to 
process  window  mode  pulse  doppler  signals.  For  this 
case,  the  digital  I  and  Q  words  are  stored.  The  phase 
and  MLM  voltage  values  are  then  computed  in  non -real¬ 
time  and  stored  on  paper  tape.  The  data  can  be  read  into 
the  MLM  registers  using  the  paper  tape  reader  of  the  hy¬ 
brid  correlator. 

LINEAR-STEP  FM  WAVEFORMS 

The  hybrid  correlator  has  also  been  used  to  pro¬ 
cess  linear-step  FM  waveforms.  In  this  mode,  the  100 
pulses  are  transmitted  at  different  frequencies.  The  fre¬ 
quency  of  each  successive  pulse  is  a  fixed  frequency  in¬ 
crement  greater  than  the  previous  pulse.  The  receiver 
shown  in  Fig.  32  is  used  but  now  the  frequency  select 
switch  is  changed  for  each  pulse.  For  a  stationary  point 
target  in  a  noise-free  environment,  a  linear  phase  slope 
across  the  return  echoes  exist  that  is  a  function  of  the 
target's  position  within  the  range  cell.  When  this  phase 
slope  is  impressed  across  the  MLM,  the  output  lobe  posi¬ 
tion  indicates  the  target's  fine  range.  For  a  complex  sta¬ 
tionary  target  in  a  real-life  environment,  the  output  is  the 
correlation  of  all  phase  slopes  (fine  ranges)  with  the  input 
data.  The  resulting  signature  can  be  used  to  identify  the 
target. 


When  this  linear-step  FM  waveform  is  used  for  a 
moving  target,  the  return  phase  is  a  function  of  both  the 
target  velocity  and  fine  range.  Since  the  doppler  fre¬ 
quency  is  a  function  of  the  transmit  frequency  (fd  = 
Vft/2C),  a  constant-velocity  point  target  will  not  generate 
a  linear  phase  slope  for  a  wide -bandwidth  transmission. 
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The  distortion  can  be  compensated  byiappropriately  ad¬ 
justing  the  time  of  transmission  of  each  ;;  Jse.  For  a 
linear-step  FM  waveform,  the  frequency  of  the  xfi1  pulse 
is  given  by: 

fn  =  tfQ  +  (n  -  1)  Af],  (28) 

where  n  =1,  2,  . ..,  100, 

fQ  =  the  transmit  frequency  of  the  first  pulse,  and 
Af  =  the  frequency  increment  between  pulses. 


To  compensate  for  the  wide  bandwidth  doppler  distortion, 
the  transmission  time  for  the  n**1  pulse  is  given  by: 


t  = 
n 


(  lo 

If  +  (n  -  l)  Af 
'  o 

\  o 


nAt 


nAt  , 


(29) 


where  At  =  the  pulse  spacing  at  the  frequency  fQ. 


REAL-TIME  SEARCH  MODE  PROCESSOR 

The  high  speed  of  the  MLM  allows  it  to  be  used  as 
a  real-time  signal  processor.  A  block  diagram  of  the 
real-time  receiver  for  the  hybrid  correlator  is  shown  in 
Fig.  33.  Associated  with  each  row  of  the  MLM  is  a  sepa¬ 
rate  phase  detector,  A/D  converter,  shift  register  delay 
line,  D/A  converter,  and  amplifier.  The  shift  register 
delay  lines  for  each  channel  are  used  for  signal  storage 
and  registration  so  that  at  any  given  time  the  phase  values 
across  the  MLM  rows  correspond  to  the  phase  values  for 
one  range  cell.  A  simple  change  in  the  clock  frequency 
allows  the  system  to  be  operated  at  data  rates  from  DC 
up  to  the  frequency  at  which  the  modulator  response  time 
becomes  a  limiting  factor. 
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Since  the  MLM  can  accept  a  new  set  of  input  data 
approximately  every  microsecond,  a  new  set  of  output 
data  is  available  every  microsecond.  A  high-speed  linear 
diode  array  is  necessary  in  the  putput  plane  in  order  to 
read  the  output  at  an  equivalent  rate.  Associated  with 
each  diode  would  be  a  threshold  detector  and  buffer 
register. 


TO  MLM  ROWS 


Fig.  33  BLOCK  DIAGRAM  OF  A  REAL-TIME  RECEIVER  FOR  THE  HYBRID  CORRELATOR 
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9.  TWO-DIMENSIONAL  RADAR  SIGNAL  PROCESSING 


Just  as  the  optical  system  can  be  used  to  perform 
two-dimensional  beam  steering  (Section  6),  the  hybrid 
correlator  can  also  be  used  to  perform  two-dimensional 
correlation  of  radar  signals.  Consider  a  frequency  di¬ 
versity  pulse  doppler  waveform  consisting  of  100  pulses 
where  the  frequency  of  the  m*h  pulse  is  given  by: 

fm  =  f0  +  n(m)  AS  .  (30) 

In  this  case  the  frequency  ordering  is  no  longer  linear  (as 
was  the  case  in  linear-step  FM  discussed  in  Section  8). 
However,  each  frequency  in  the  set  of  100  pulses  is  dif¬ 
ferent,  so  that  the  whole  bandwidth  is  covered  uniformly. 
For  a  point  target  in  a  noise-free  environment,  a  linear 
phase  slope  with  respect  to  frequency  exists  across  the 
return  echoes  that  is  a  function  of  the  target's  fine  range. 
Also  a  linear  phase  slope  with  respect  to  time  exists  across 
the  return  echoes  that  is  a  function  of  the  target's  velocity 
(ignoring,  for  now,  the  doppler  frequency  dependence  on 
transmit  frequency).  In  other  words,  a  two-dimensional 
phase  slope  exists  across  the  frequency /time  plane.  The 
function  of  a  correlator  is  to  determine  this  two-dimen¬ 
sional  phase  slope  and  therefore  the  velocity  and  fine  range 
of  the  target.  The  two-dimensional  phase  slope  corre¬ 
sponding  to  a  given  velocity  and  fine  range  is  assumed. 

The  correlation  of  the  return  echoes  with  this  phase  slope 
is  given  by: 


M 

C(V,  R)  =  Z  A  exp  (to)  exp(-  exp( n(m)  RH  (31) 

m  =  l  At 

where  =  the  amplitude  of  the  mth  echo, 

<p  -  the  phase  of  the  m^  echo. 
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V  =  the  assumed  velocity  of  the  target. 

At  =  the  time  increment  between  pulses, 

Xt  =  the  wavelength  of  the  transmitted  waveform-, 

R  =  the  assumed  fine  rahge  of  the  target, 

Af  -  the  frequency  increment,  and 
c  =  the  velocity  of  light.  i  : 

Notice  that  this  equation  is  similar  in  form  to  Eq.  (24), 
which  describes  the  transform  pattern  of  a  thinned  ahray. 
When  the  associations  are  made: 

AX  -  At  , 

e  -  v.  ,  .  '  • 

AY  -  Af  ,  and  .  1  . 

7}  R  ,  (32) 

Eq.  (24)  can  be  interpreted  as  a  two-dimensional  corre¬ 
lation  function,  the  phase  of  the  mth  echo  is  impressed 
on  the  mth  radiator-,  so. that  the  spacing  of  the  radiators  • 
in  the  X  direction  corresponds  to  the  relative  transmit 
time  of  the  pulses  and  the  spacing  of  the  radiators  in  the 
Y  direction  corresponds  to  the  relative  frequency  o.f  the 
pulses.  The  position  along  the  £  axis  in  the  transform 
plane  corresponds  to  the  velocity  of  the  target,  and  the 
position  along  the  r\  axis  in  the  transform  plane  corre¬ 
sponds  to  the  fine  range  of  the  target.  The  optical  sys¬ 
tem  not  only  performs  the  correlation  instantaneously, 
but  also  the  correlation  is  performed  simultaneously  for 
all  velocities  and  fine  ranges.  Figure  34  depicts  this 
two-dimensional  signal  processing  capability  of  the  hy¬ 
brid  correlator. 

The  doppler  frequency  dependence  on  transmit : 
frequency  can  again  be  compensated  for  by  appropriately 
altering  the  transmit  time  of  the  pulses.  The  new  trans¬ 
mission  time  for  the  mth  pulse  is  given  by: 
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t  =  m  ( 7 — —7 — r-r  )  At  ,  (33) 

m  \  f  +  n(m)Af  /  o 
'  o 

where  At  =  the  time  increment  at  f  . 

o  o 

When  the  time -frequency  code  of  the  radar  wave¬ 
form  is  determined  judiciously,  the  velocity  and  fine 
range  resolution  of  the  full  time -bandwidth  can  be  main¬ 
tained.  This  is  analogous  to  the  selection  of  an  appro¬ 
priate  mask  for  two-dimensional  beam  steering  which 
maintains  the  resolution  of  the  entire  array.  Also,  the 
sidelobes  that  are  obtained  from  the  two-dimensional 
beam-steering  masks  are  the  sidelobes  associated  with 
the  ambiguity  diagram  of  the  radar  waveform.  Thus, 
the  ’'best"  masks  that  were  obtained  for  two-dimensional 
beam  steering  will  also  be  the  best  time -frequency  code 
for  the  radar  waveforms.  The  successful  two-dimensional 
beam-steering  results  (Section  6)  indicate  the  feasibility  of 
using  the  hybrid  correlator  for  two-dimensional  radar  sig¬ 
nal  processing. 

For  window-mode  processing,  the  receiver  shown 
in  Fig.  32  is  used  again.  Now,  the  frequency  selector 
is  switched  to  the  appropriate  frequency  for  each  pulse. 

The  ultimate  capability  of  the  hybrid  correlator 
would  be  to  process  two-dimensional  waveforms  in  a 
real-time  search  mode.  Since  the  MLM  can  accept  a 
new  set  of  input  data  every  microsecond,  a  new  set  of 
correlated  data  occurs  every  microsecond.  The  corre¬ 
lator  is  now  limited  by  the  speed  at  which  the  two- 
dimensional  field  of  output  data  can  be  interrogated.  A 
suitable  means  is  not  currently  available  for  interrogat¬ 
ing  the  two-dimensional  data  plane  at  megahertz  rates. 

A  high-speed  two-dimensional  array  of  photodiodes 
followed  by  threshold  circuits  and  amplifiers  appears 
to  offer  the  most  promise. 
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APPENDIX  A 

DEFLECTION  CHARACTERISTICS  OF  THE  MEMBRANE 
LIGHT  MODULATOR 


SUMMARY 

A  mathematical  model  for  the  deflection  character¬ 
istics  of  the  Membrane  Light  Modulator  is  described  in  this 
appendix.  A  power  series  solution  has  been  obtained  which 
describes  the  static  deflections  of  the  membrane  as  a  func¬ 
tion  of  applied  voltage  for  small  deflections.  The  signifi¬ 
cant  conclusions  obtained  from  this  analysis  are: 

1.  Deflection  is  proportional  to  the  square  of  the 
applied  voltage  {for  small  deflections  only). 

2.  Maximum  stable  deflection  is  44%  of  the  initial 
separation  and  is  independent  of  the  dimen¬ 
sions. 

3.  If  it  is  desired  to  utilize  a  25%  deflection,  a 
further  increase  of  only  10%  in  voltage  will 
cause  instability  {beyond  which  the  membrane 
will  pull  down  onto  the  fixed  electrode). 


MATHEMATICAL  MODEL 

The  cross-section  of  an  MLM  element  is  shown  in 
Fig.  A-l.  The  collodion  membrane  which  is  coated  with 
a  conducting-reflecting  coating  is  stretched  over  a  cylin¬ 
drical  well  of  radius  R0  and  depth  D.  The  initial  tension 
of  the  membrane  is  T.  The  element  possesses  circular 
symmetry  about  the  z  axis  and  therefore  cylindrical 
coordinates  are  used. 

When  a  voltage  V  is  applied  between  the  mem¬ 
brane  and  the  electrode,  the  membrane  deflects  to  a  new 
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ILLUMINATION 


ANTIMONY  COATE 0 


ELECTRODE 


R0 


Fig.  A-1  CROSS  SECTION  OF  AN  MLM  ELEMENT 
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static  position.  The  maximum  deflection  (A)  occurs  at  the 
center  of  the  element.  We  are  interested  in  finding  the 
deflection  A  as  a  function  of  the  applied  voltage  V. 

Four  assumptions  are  made: 

1.  The  membrane  has  negligible  stiffness;  i.  e. , 
only  tension  acting  tangentially  need  be  con¬ 
sidered. 

2.  For  the  small  D/R0  ratios  in  the  practical 
case,  the  tension  T  may  be  considered  a 
constant. 

3.  For  these  small  D/R0  ratios,  dz/dr  is  small 
compared  to  l,  permitting  a  simplified  ex¬ 
pression  for  the  curvature. 

4.  Electrostatic  fringing  or  end  effects  may  be 
neglected. 

From  the  "hoop-tension"  formula  in  mechanics,  a 
membrane  is  in  equilibrium  when  the  product  of  tension 
and  curvature  is  equal  to  the  radial  pressure.  The  radial 
pressure  in  this  case  is  given  to  a  good  approximation  by 
the  parallel  plate  capacitor  expression,  since  electro¬ 
static  fluxlines  are  always  normal  to  the  electrodes. 
Therefore,  at  any  point  on  the  membrane,  the  electro¬ 
static  pressure  is  (in  MKS  practical  units): 

v2 

p  =  — newtons  per  square  meter,  (A  - 1 ) 

2z2 

where  V  =  voltage  applied  between  the  membrane  and  the 
bottom  electrode 

z  =  electrode  separation  in  meters 

1  -9 

€  =  permittivity  ~  x  10  F/m. 
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The  membrane  possesses  a  double  curvature.  The 
total  curvature  in  such  a  case  is  the  sum  of  the  curvature 
in  orthogonal  directions.  In  rectangular  coordinates, 
curvature  along  the  x  direction  is  given  by:_ 

curvature  = - (A -2) 

1  +  <dz/dx) 


For  a  relatively  flat  surface,  the  denominator  is  approxi¬ 
mately  unity  so  that  d2z/dx2  suffices.  Thus,  the  total 
curvature  for  a  surface  that  is  curved  with  recoct  to  y 
as  well  as  x  is  simply: 


total  curvature  = 


(A-3) 


The  corresponding  expression  in  cylindrical  coordinates 
may  be  determined  by  noting  that  Eq.  (A-3)  is  the 
Laplacian  operator  or  V2  in  rectangular  coordinates.  Its 
equivalent  in  cylindrical  coordinates  is: 


curvature  = 


I JL 

r  dr 


(r 


(A-4) 


In  our  model,  due  to  symmetry,  there  is  no  variation  with 
0.  Hence, 


curvature  = 


Id.  dzv 
r  dr  (r  dr 


1  dz 
r  dr 


(A-5) 


For  a  membrane  tension  of  T  N/m,  the  Mhoop  ten¬ 
sion"  equation  becomes: 


d2z  1  dz  _  V2t  _JK 

,  2  r  dr  “  _  2  '  2 

dr  2Tz  z 


where  K  =  V“c/  2T. 


(A-6) 
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SOLUTION 


The  solution  of  Eq.  (A -6)  will  give  the  membrane  de¬ 
flection  versus  r  for  various  values  of  K.  This  is  a  non¬ 
linear  differential  equation  with  no  obvious  closed  form 
solution.  The  first  five  terms  of  the  power  series  solution 
of  Eq.  (A-6)  is: 


..]  .  (A-  7) 


The  coefficient  of  the  last  term  is  approximate.  There  is 
no  recognizable  closed  form  which  will  give  this  succession 
of  coefficients,  an  expected  situation  in  view  of  the  com¬ 
plexity  of  the  differential  equation.  When  the  boundary 
condition  that  z  =  z0  at  r  =  0  is  inserted  into  Eq.  (A -7),  the 
arbitrary  coefficient  is  found  to  be  a  -  z0  (where  z0  =  D  -  A). 


What  we  are  particularly  interested  in  is  not  the 
equation  for  the  full  membrane,  but  rather  the  spacing  at 
the  center  (zQ).  We  may  therefore,  insert  into  Eq.  (A -7) 
the  boundary  condition  z  =  Dat  r  =  Rq: 


KR‘ 


D  =  z  [1  + 
o  1 


4z 


_l_ 

32 


1 

520 


(A-8) 


This,  then  is  the  final  equation  relating  zQ  to  K  for  a  given 
R0  and  D. 


NUMERICAL  RESULTS 

From  Eq.  (A-8),  it  is  seen  that  the  series  (within 
the  brackets)  is  fully  determined  by  the  value  of  KRq/zq3. 
Hence,  one  approach  to  numerical  results  is  to  assume  a 
value  for  KRs/zQ3,  sum  the  terms  in  the  series,  and  then 
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solve  for  zQ.  From  the  assumed  value  of  KRq/z0^  and  the  , 
solved  value  of  zD,  K  may  now  be  determined.  This  pro¬ 
cess  works  well  for  small  values  of  KR§/z03  (correspond¬ 
ing  to  a  large  zQ  or  small  deflection)  where  the  series 
converges  rapidly.  However,  for  larger  deflections,  five 
terms  in  the  series  are  insufficient1  for  accurate  results 
and  the  method  fails.  Even  if  one  were  to  expend  the  con¬ 
siderable  labor  necessary  to  obtain  a  few  additional  coeffi¬ 
cients,  it  would  still  be  impossible  to  determine  the  maxi¬ 
mum  stable  deflection  with  any  precision.  Consequently, 
a  computer  numerical  solution  of  Eq.  (A -6)  was  sought, 
using  the  boundary  condition  D  =  1  j*m  and  R0  =  18.  5  jim. 

The  computer  approach  was  as  follows:  A  value  of 
K  is  assumed,  and  a  corresponding  value  for  z0  is  guessed. 
From  symmetry,  the  first  derivative  of  z  at  r  =  o  is 
z0'  =  o.  From  Eq.  (A -7),  the  second  ^derivative  is 
z0"  =K /  2zQ2.  Using  these  initial  values,  the  differential 
equation  is  solved  numerically  to  give  a  value  for  z  at 
r  =  R  =  18.  5  jim.  If  this  value  is  other  than  z  =  1  jxm,  a 
new  value  for  zQ  is  guessed.  Finally  a  value  for  zQ  is  ob¬ 
tained  for  which  z  =  1  at  r  =  18.  5.  This  process  is  re¬ 
peated  for  increasing  values  of  K  until  a  value  for  K  is 
reached  where  the  boundary  condition  is  satisfied  over  an 
interval  of  zQ.  At  this  point,  the  curve  of  deflection  ver¬ 
sus  K  is  vertical  and  is,  therefore,  the  limit  of  stability. 
For  attempted  solutions  for  greater  values  of  K,  it  is 
found  that  z  >  1  at  r  =  18.  5  regardless  of  the  choice  of  zD. 

A  plot  of  deflection  (A)  is  shown  in  Fig.  A -2. 

Equation  (A -8)  may  be  rewritten  in  the  following 

form: 
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It  is  evident  that  the  ratio  z0/D  (and,  therefore,  the  frac¬ 
tional  deflection  A/D)  is  completely  determined  by  KR^/D3. 
Therefore,  the  plot  of  Fig.  A -2  i3  valid  for  other  values  of 
R0  and  D  (provided  the  assumptions  are  still  satisfied) 
if  the  ordinate  is  now  considered  to  be  A/D  and  the  abscissa 
K(R/18.  5)2/D3. 
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APPENDIX  B 

HYBRID  CORRELATOR  ELECTRONICS 


A  block  diagram  of  the  Hybrid  Correlator  elec¬ 
tronics  has  been  described  in  Section  4  (Fig.  12).  A  more 
detailed  description  of  some  of  the  components  in  this 
block  diagram  will  be  given  in  this  appendix  including  the 
protection  electronics  that  were  fabricated  to  keep  the 
MLM  elements  from  bottoming  out.  The  modifications  of 
the  MLM  electronics  required  for  beam  steering  with  the 
crystal  will  also  be  discussed.  Finally,  a  description  of 
the  scanning  mirror  readout  system  will  be  given. 


MLM  BEAM-STEERING  ELECTRONICS 

A  linear  phase  slope  must  be  impressed  across  the 
MLM  rows  in  order  to  steer  the  beam  of  this  phased  array 
antenna  in  one  dimension.  The  Beam  Position  Generator 
(BPG)  or  a  digital  computer  can  be  used  to  generate  the 
appropriate  voltages  for  each  row  of  the  MLM  in  order  to 
obtain  these  linear  phase  slopes.  The  MLM  phase  depen¬ 
dence  on  the  square  of  the  applied  voltage  must  be  con¬ 
sidered  in  making  these  calculations. 

An  M-460  digital  computer  can  be  used  to  generate 
the  voltages  required  to  steer  the  beam  into  the  100  Ray¬ 
leigh  resolvable  beam  positions  of  the  MLM  array.  The 
voltage  words  are  loaded  sequentially  into  the  input  regis¬ 
ters  of  the  correlator  electronics  in  the  following  manner. 
First  a  data  request  signal  is  sent  to  the  computer  (via 
the  computer  control  electronics  in  Fig.  12).  When  ready, 
the  computer  returns  a  data  acknowledge  pulse  to  the  com¬ 
puter  control  box  and  data  for  the  first  row  of  the  MLM  to 
the  input  register.  The  data  acknov/ledge  pulse  is  used  to 
strobe  the  data  into  the  registers.  After  the  100th  data  re¬ 
quest  and  data  acknowledge,  an  "external  function"  pulse 
is  sent  by  the  computer  which  is  used  to  transfer  the  data 
into  the  storage  register  and  clear  the  input  register. 
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An  interrupt  line,  sn  interrupt  acknowledge  line, 
and  a  number  of  data  lines  are  also  connected  to  the  com¬ 
puter  for  the  purpose  of  communicating  with  the  computer. 
The  data  lines  are  set  to  a  desired  code  (via  front  panel 
toggle  switches)  which  will  tell  the  computer  the  desired 
transmission  mode. 

As  discussed  in  Section  8,  the  digital  computer  is 
used  off-line  to  determine  the  phases  (and  also  the  re¬ 
quired  MLM  voltages)  of  the  radar  pulses.  Radar  data 
from  the  computer  can  be  accepted  continuously.  With  a 
computer  output  rate  of  one  word  every  50  jxs.  it  will  take 
5  ms  to  load  the  MLM  for  one  set  of  data.  A  "stop"  switch 
is  also  available  which  allows  a  set  of  data  to  be  held  on 
the  MLM  for  detailed  study. 

The  data  generated  by  the  M-460  can  be  recorded 
on  7-level  paper  tape.  The  first  6  bits  are  used  for  the 
MLM  voltage  data  and  bit  7  is  used  for  the  load  pulse  (ex¬ 
ternal  function  of  computer).  100  lines  of  data  on  the  paper 
tape  are  required  corresponding  io  the  100  rows  of  the  MLM. 
Bit  7  is  punched  in  line  101.  This  bit  is  used  to  load  the 
data  into  the  output  register  (ano  clear  the  input  register) 
nnd  also  can  be  used  to  stop  the  tape  reader  after  each 
block  ot'  data.  The  other  6  bits  in  this  line  are  used  as  a 
code  for  the  next  set  of  data. 

The  control  logic  for  the  paper  tape  reader  was 
built  to  allow  a  great  deal  of  flexibility  in  the  mode  of 
operation.  The  modes  that  arc  used  most  often  are  listed 
below: 


1.  Tape  Run.  In  the  tape  run  mode,  the  reader 
runs  continuously  at  its  maximum  speed  of  625  lines  per 
second  (which  results  in  approximately  6  blocks  of  data 
per  second).  It  can  be  made  to  stop  after  any  of  the  blocks 
of  data  or  at  any  of  the  codes.  (The  desired  code  is  se¬ 
lected  via  toggle  switches  on  the  front  panel. ) 
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2.  Block  Rate  Control.  The  rate  at  which  blocks 
are  read  by  the  tape  reader  can  be  controlled  in  the  Block 
Rate  Control  Mode.  The  block  rate  can  be  varied  from  6 
blocks/second  to  1  block  every  10  seconds.  In  this  mode 
the  100  lines  of  a  block  are  read  at  a  625  line  per  second 
rate  and  then  the  reader  stops  at  the  end  of  the  block  for 
the  desired  length  of  time  before  reading  the  next  block  of 
data.  This  block  reading  can  be  stopped  after  any  block 
or  at  any  of  the  code  words. 

3.  Single  Block  Load.  The  "momentary"  switch 
allows  single  blocks  to  be  loaded. 

As  discussed  in  Section  6,  a  PL/I  computer  pro¬ 
gram  has  been  written  (for  use  on  the  IBM  360/91)  to  gen¬ 
erate  the  voltages  required  for  steering  the  beam  from 
the  thinned  MLM  arrays  in  two-dimensions.  The  output 
from  this  computer  program  is  punched  on  paper  tape  in 
a  format  which  is  compatible  with  the  paper  tape  reader. 

In  order  to  provide  a  fast  beam-stepping  mode,  a 
special  purpose  digital  package,  called  the  Beam  Position 
Generator  (BPG)  in  Fig.  12,  has  been  fabricated  which  can 
be  used  to  load  data  into  the  input  registers  in  a  quasi¬ 
parallel  mode.  For  simplicity  in  the  digital  design  of  the 
BPG,  the  MLM  is  considered  to  have  128  resolvable  beam 
positions.  These  128  positions  are  defined  as  the  unit 
beam  positions  and  are  further  divided  into  fractional 
positions  of  1/2,  1/4,  and  1/8.  As  shown  in  Fig.  B-l,  the 
phase  of  all  of  the  rows  can  be  determined  by  adding  the 
appropriate  offset  to  the  phases  of  the  first  8  rows.  The 
phase  increment  (A<p)  is  a  function  of  the  desired  beam 
position.  Using  this  technique  to  generate  the  phase  slope, 
the  input  register  is  loaded  in  13  parallel  groups.  Each 
group  contains  8  voltage  words  (except  for  the  last  group 
which  contains  4)  which  are  loaded  sequentially  into  the 
input  registers.  In  this  mode  only  8  clock  pulses  are  re¬ 
quired  to  load  the  input  registers  instead  of  the  100  clock 
pulses  required  in  a  purely  sequential  operation.  The 
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stepping  time  from  one  beam  position  to  the  next  is  limited 
by  the  digital  logic  to  around  2  ^s.  The  price  for  this 
speed  capability  is  the  requirement  for  13  square  root 
coders  and  13  sets  of  output  data  from  the  BPG. 

The  BPG  is  useful  for  evaluating  the  MLM  system 
without  the  need  of  the  computer.  It  has  several  modes  of 
operation: 

1.  Free-Running.  The  beam  from  the  MLM  is 
stepped  continuously  through  all  the  beam  positions  in 
full  or  fractional  steps.  The  stepping  speed  is  controlled 
by  means  of  the  basic  clock  frequency. 

2.  Selected  Beam  Position.  The  simulator  runs 
until  it  reaches  a  selected  beam  position  and  then  stops. 
The  desired  beam  position  number  is  selected  by  means 
of  toggle  switches. 

3.  Toggle  Mode.  This  is  a  high  speed  mode  used 
to  test  the  frequency  response  of  the  MLM  which  by¬ 
passes  the  2  per  step  limitation  in  the  free-running 
mode.  The  selected  beam  position  is  stored  in  the  input 
register.  The  output  register  then  rapidly  alternates  be¬ 
tween  this  selected  beam  position  and  zero  at  each  clock 
pulse. 

4.  Noise  and  Offset  Mode.  This  mode  allows 
digital  noise  (from  an  external  source)  or  a  uniform  phase 
offset  to  be  added  sequentially  to  the  simulator-generated 
phase  words.  In  order  to  obtain  noise  that  is  independent 
row  to  row,  the  levels  for  each  row  must  be  determined 
sequentially.  This  requires  100  clock  pulses  which  re¬ 
sults  in  a  low-speed  operation. 


MLM  PROTECTION  CIRCUITS 

As  described  in  Appendix  A,  when  overvoltages  are 
applied  to  the  modulator  elements,  they  bottom  out.  The 
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collodion  membrane  adheres  to  the  SiO  protective  layer 
over  ‘.he  electrode  structure.  A  voltage  clamp  cir  cuit  has 
been  fabricated  which  prohibits  the  outputs  from  any  of 
the  MLM  amplifiers  from  exceeding  a  certain  threshold 
value  (usually  set  just  above  the  2tr  phase  voltage).  Any 
excess  current  due  to  an  amplifier  failure  or  a  turn-on 
or  turn-off  transient  is  drained  off  by  this  clamp  circuit. 

The  modulator  elements  were  also  observed  to 
bottom  when  exposed  to  relatively  low  voltages  for  a  long 
period  of  time.  A  special  circuit  has  been  built  to  pre¬ 
vent  a  voltage  from  staying  on  an  element  for  longer  than 
10  seconds.  If  a  load  pulse  does  not  occur  at  the  output 
register  every  10  seconds  or  less,  this  circuit  automati¬ 
cally  resets  the  output  registers  to  a  zero  level.  This 
circuit  works  independent  of  the  mode  of  operation  of  the 
electronics  (i.  e. ,  computer,  paper  tape,  etc.).  A  switch 
allows  this  circuit  to  be  disconnected  during  the  optical 
alignment  procedure. 


CRYSTAL  BEAM-STEERING  ELECTRONICS 

The  paper  tape  mode  of  the  MLM  electronics  is 
used  to  beam  steer  with  the  crystal.  As  discussed  in 
Section  5,  a  PL/ 1  computer  program  has  been  written 
(for  use  on  the  IBM  360/91)  which  generates  the  voltages 
required  to  steer  the  beam  of  the  crystal  array.  This 
program  includes  the  offset  phases  required  to  compen¬ 
sate  for  the  phase  distortion  across  the  crystal  surface. 
The  output  of  this  program  is  punched  on  paper  tape  in  a 
format  which  is  compatible  with  the  paper  tape  reader  of 
the  hybrid  correlator.  Since  there  are  only  46  crystal 
elements,  only  the  first  46  input  and  output  registers  of 
the  hybrid  correlator  are  used. 

A  potentiometer  is  used  at  the  output  of  each  of 
the  MLM  amplifiers  in  order  to  compensate  for  the  dif¬ 
ferent  2 it  phase  voltages  of  the  MLM  (83  volts)  and  the 
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crystal  (around  32  volts).  The  phase  modulation  sensitivity 
of  the  individual  crystal  elements  is  not  exactly  the  same 
and  therefore  these  potentiometers  can  also  be  used  to  ad¬ 
just  the  outputs  to  the  different  2ir  voltages  of  each  element 
(which  range  from  30  to  33  volts). 


SCANNING  MIRROR  READOUT  SYSTEM 

The  readout  system  used  to  read  the  data  from  the 
output  plane  of  the  hvbrid  correlator  is  shown  in  Fig.  B-2. 
The  scanning  mirror  sweeps  the  area  of  interest  in  the 
readout  plane  across  the  face  of  a  photomultiplier  tube. 
Thus  the  current  amplitude  output  of  the  photomultiplier 
as  a  function  of  time  corresponds  to  the  intensity  of  the 
light  in  the  readout  plane  as  a  function  of  distance  along 
the  scan  interval.  The  output  from  the  photomultiplier  is 
displayed  on  an  oscilloscope. 

A  projection  lens  is  used  to  focus  the  readout  data 
plane  via  the  scanning  and  fixed  mirrors  onto  the  slit  in 
front  of  the  photomultiplier  tube.  Due  to  the  small  size 
of  the  scanning  mirror  (i  by  i  inch)  it  was  necessary  to 
use  a  short  focal  length  projection  lens.  Since  the  image 
is  many  focal  lengths  away  from  the  lens,  the  readout 
plane  is  approximately  in  the  back  focal  plane  of  the  pro¬ 
jection  lens.  Therefore,  in  the  front  focal  plane  is  the 
Fourier  transform  of  the  readout  data  plane  (which  is 
just  a  scaled  image  of  the  masked  modulator).  The  size 
of  this  image  is  approximately  the  ratio  cf  the  focal 
lengths  of  the  projection  lens  and  the  transform  lens 
(/2//1)*  The  scanning  mirror  is  placed  in  this  plane. 

Since  :he  image  is  smaller  than  the  size  of  the  mirror, 
all  the  orders  that  are  passed  through  the  projection  lens 
can  be  scanned  past  the  photomultiplier.  The  fixed  mirror 
is  used  only  as  a  matter  of  convenience  to  allow  the  photo¬ 
multiplier  to  be  mounted  in  an  uncluttered  location. 

The  scanner  utilizes  the  model  G-125  galvanometer 
scanning  mirror  driven  by  the  AX-100  driver  amplifier 
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Fig.  B-2  SCANNING  MIRROR  READOUT  SYSTEM 
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manufactured  by  General  Scanning,  Inc.  The  galvanometer 
mirror  can  scan  up  to  10®  maximum  peak-to-peak  deflection 
and  has  a  resonant  frequency  of  approximately  1*400  Hz.  A 
linear  scan  was  obtained  by  using  the  sawtooth  output  from 
a  Tektronix  oscilloscope  as  a  signal  source  for  the  driver 
amplifier.  Scan  frequencies  as  high  as  100  Hz  have  been 
used. 
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APPENDIX  C 

ELECTRO-OPTIC  PROPERTIES  OF  THE  LITHIUM 
TANTALATE  MODULATOR 


It  has  been  shown  by  many  authors  (Refs  3,  7, 
and  8)  that  the  relative  optical  impermeability  tensor 
(and  therefore  the  refractive  index)  of  certain  crystals  is 
a  function  of  the  applied  electric  field.  The  impermea¬ 
bility  tensor  0^  f°r  a  clamp  crystal  (clamping  is  neces¬ 
sary  in  order  to  eliminate  the  effect  due  to  piezoelec¬ 
tricity  and  photoelasticity)  can  be  described  by  the  tensor 
equation; 


0..  =  0.,  +  r  .  E.  +  p..  .  E.  E 
lk  lk  ik]  3  lkjm  3  m 

(i,  3,  k,  m  =  1,  2,  3) 


(C-l) 


where  /3.k  =  tensor  describing  the  relative  impermea¬ 

bility  for  zero  field, 

E.  =  tensor  describing  the  applied  electric 
3  field  (E), 

r  =  tensor  of  rank  three  whose  27  components 
1  J  are  the  linear  electro -optic  coefficients 

("Pockels"  effect), 

p.,  .  -  tensor  of  rank  four  whose  components  are 

J  the  quadratic  electro-optic  coefficients 

("Kerr"  effect). 

By  convention,  X3  is  chosen  parallel  to  the  ac 
axis  (optic  axis),  Xj  is  parallel  to  the  a  axis^and  X2  is 
parallel  to  the  6  axis. 

For  small  applied  fields  (such  as  those  used  in 
this  experiment),  the  Kerr  effect  is  negligible  with  respect 
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to  the  Pockets  effect.  The  27  components  of  the  ten¬ 
sor  can  be  reduced  to  18  by  using  symmetry  arguments. 
The  resulting  equation  is  customarily  written  in  the  fol¬ 
lowing  manner: 


~*n  ■  < 

r  r  r 

11  12  13 

022  -  0°2 

r  r  r 

21  22  23 

033  -  03 

r  r  r 

31  32  33 

X 

^23 

r4l  r42  r43 

E2 

^31 

r51  r52  r53 

A 

_  _ 

r61  r62  r63_ 

where  the  electro -optic  coefficients  r^  are  now  repre¬ 
sented  by  a  3  x  6  matrix. 


Since  lithium  tantalate  is  a  trigonal  crystal  belong¬ 
ing  to  the  3m  point  group,  it  is  optically  uniaxial  (when  no 
field  is  present)  and  therefore  0^  =  02.  Also,  the  r^ 
matrix  reduces  to  (Ref.  5): 


ik 


0 

0 

0 

0 

r 

-r 


51 

22 


“r22  ri3 


r22  ri3 
0  r 


33 


51 


(C-3) 


where  r22  -  -r12  =  -r6l;  r51  =  r42;  and  r^  =  r23. 
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The  impermeability  matrix  is1  found  by  substituting 
Eq.  (03)  into  Eq.  (02): 


^1  ‘  r22E2  +  ri3E3  . 


-r  E 
22  1 


r  E 
51  1 


~r22El 


^l+r22E2  +  ri3E3  r51E2 


r  E 
51  1 


r  E 
51  2 


^3  *  r33E3 


(04) 


a 

(U)  When  the  electric  field  is  applied  parallel  to  the  X3 
axis  (Ej  *  E2  =  0),  the  impermeability  matrix  reduces  to: 


*l+ri3E3 


*1  +  r!3E3  ■  ° 


^3  +  r33E3 


(C— 5) 


The  8  ellipsoid  (Q..  X.  X.  «  1)  is  described  by: 
i  1  3  1 


K+r!3E3)Xl2  +  ^l+ri3E3)X22 


(C-6) 


+  <*3  +  r33E3>  X32  =  1  • 


Notice  that  this  is  still  in  principle  axis  form  which 
means  that  a  field  applied  parallel  to  the  X3  axis  does  not 
rotate  the  principle  axes.  The  effect  of  the  field  is  to 
change  the  length  of  the  axes.  Since  the  coefficients  of  the 
Xj2  and  X22  terms  are  equal  (£°  +  r13E3),  the  crystal  re¬ 
mains  uniaxial. 


For  light  polarized  in  the  X3  direction  and  propa¬ 
gating  in  the  X2  direction,  the  phase  modulation  (Atp)  rela¬ 
tive  to  a  crystal  of  the  same  length  but  with  zero  applied 
field  is: 
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3**3  =  -U  3An3  , 


(C-7) 


where  1 2  =  length  of  crystal  in  the  Xg  direction, 

A 

3  3  =  change  of  index  of  refraction  in  the  Xg  di¬ 

rection  due  to  the  applied  field  (E3) 

=  n3  - 

X  =  wavelength  of  the  illuminating  beam. 

The  impermeability  of  a  crystal  is  related  to  the 
index  of  refraction  by 


0  =  1/n  . 


(C-8) 


Therefore: 


3^  3  =  3^3  -  =  r33E3 


1  2  2 
(n  y  (n  S 


(C-9) 


Rearranging  Eq.  (C-9)  yields: 


r33E3  = 


<n3  "n3)(n3  +  n3) 

.  \2  .  \2 

(n3)  (n3) 


(C-10) 


For  —  «  1,  n3  ftsng  and  Eq.  (C-10)  reduces  to; 
n3 


0An0  ‘n3  r33E3 
J  .5  2 


(C-U) 
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Substituting  into  Eq.  (07)  and  letting  E3  =  V^/lg, 
the  phase  modulation  due  to  ar  applied  field  parallel  to  the 
Xg  direction  of  light  polarized  parallel  to  the  Xg  direction 
and  propagating  in  the  X2  direction  is: 


l2  n3  r33  V3 
XI. 


(012) 


Similarly,  for  light  polarized  parallel  to  the  X^  direction 
and  propagating  in  the  X2  direction,  the  phase  modulation 
due  to  an  applied  field  parallel  to  the  X3  direction  can  be 
shown  to  be: 


-it  t2  ^  r  13 
XI. 


(C- 13) 


The  Xg  and  Xj  directions  represent  the  allowed  di¬ 
rections  of  polarization  in  the  crystal  for  light  propagating 
in  the  X2  direction.  When  light  is  polarized  at  an  angle 
with  respect  to  these  axes,  the  light  (as  far  as  the  crystal 
is  concerned)  breaks  up  into  Xg  and  Xj  components.  The 
phase  modulation  through  the  crystal  for  toe  Xg  compo¬ 
nent  is  given  by  Eq.  (C-12),  and  for  the  X^  component  by 
Eq.  ( C- 13).  The  relative  phase  difference  between  these 
components  at  the  exit  surface  of  the  crystal  is  given  by: 


3^3  -  1A<P3 


_7r  i2  V3  /  „3  _3 

Xl3  \r33  3  _ri3ni 


(014) 


This,  in  general,  results  in  an  elliptically  polar¬ 
ized  beam. 
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APPENDIX  D 

OPTICAL  PHASED  ARRAY  TRANSFORM  PATTERNS 


The  equations  describing  the  Fourier  transform 
patterns  of  optical  phased  arrays  will  be  discussed  tutorialiy 
in  this  appendix.  The  reader  is  directed  to  Ref.  9  for  a 
description  of  the  Fourier  transform  properties  of  a  lens 
and  to  Ref.  10  for  a  discussion  of  optical  diffraction  theory. 

When  an  object  is  placed  in  the  front  focal  plane  or 
(X,  Y)  space  of  a  perfect  lens,  its  spatial  Fourier  trans¬ 
form  exists  in  the  back  focal  plane  or  space  and  is 

described  by  the  following  equation  (Ref,  9): 

U<4,rj)  =  J  J  U(X,  Y)  exp{  -j  ~  (X£  +  Ytj)]  dx  dy  ,  (D-l) 

-00  -0D 

where  U(X,  Y)  -  the  complex  light  distribution  in  the 
input  plane, 

U(£,T7>  =  the  complex  light  distribution  in  the 

transform  plane, 

X  *  the  wavelength  of  the  coherent  illumina¬ 

tion,  and 

f  -  the  focal  length  of  the  lens. 

Typically  the  input  function  is  described  by  an 
amplitude  term,  A(X,  Y),  and  a  phase  term,  cxp[yp(X,  Y)], 
such  that: 

U(X,  Y)  =  A(X,  Y)  exp[jp(X.  Y)]  .  (D-2) 

Consider  first  a  linear  array  of  coherent  isotropic 
light  sources  which  are  equally  spaced  along  the  X  axis 
(as  shown  in  Fig.  8).  Equation  D-l  reduces  to  one  dimen¬ 
sion  and  can  be  written  in  the  series  form: 


Preceding  page  blank 
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uaU*^)  =  E  Am  ex[j(jpm)  exp[-j  |y  (m  -  ^~i)  AX£j  ,  (D-3) 

m  =  l  A*' 


where  M  -  the  number  of  point  radiators, 

Am  =  amPlitucle  °f  the  mth  radiator, 

O  =  the  relative  phase  of  the  ni  radiator,  and 
m 

AX  -  the  center  spacing  of  the  elements. 

When  the  radiators  are  of  the  same  amplitude  and 
phase  (i.  e. ,  Am  =  A  and  <pm  =  0),  Eq.  (D-3)  reduces  to: 


ua(£)  =  S  A  exp[-j  |p(m  -^ii)  AX£]  . 


(D-4) 


Equation  (D-4)  can  be  rearranged  into  the  form: 


M  1 

U  ({)  *  A  exp[+j£(M+l)  AX£j  t  Zm-l|,  (D-5) 

V  m=0  J 


where: 


Z  =  exp[-j0(AX£)]  . 


(D-6) 


The  sum  in  the  right-hand  brackets  of  Eq.  (D-5)  is  in  the 
form  of  a  geometric  series  whose  solution  is  given  by: 


M 

Z  Z 
m=0 


m  1-Z 

1-Z 


(D-7) 


Substituting  Eqs.  (D-6)  and  (D-7)  into  (D-5)  yields: 

,t  /„*  a  exp[j(Mi7AX£)/X/l-  e:;p[-j(Mff AX£)/X/) 
Ua(t)-A  exp(jVAX?TvT-  exp(  -  jit  AX£  /  X/l 
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Using  Euler's  equation,  the  equation  describing 
the  transform  pattern  due  to  a  one-dimensional  array  of 
isotropic  radiators  is 

VO  - A  •  (D- 

The  power  in  the  transform  plane  can  be  found 
from  Eq.  (D-9): 

Pa(V  -  |ua(«)|2  =  u a(«)ua*(«) 


a2  sin2  [(Mff  AX£ )/Xf]  _  (D.10] 

sin2  [(i7£X£)/X/] 

This  equation  describes  a  one -dimensional  array 
of  grating  lobes  which  occur  whenever  both  the  numerator 
and  denominator  of  this  equation  go  to  zero  simultaneously. 
This  occurs  whenever: 


7TAX£ 

-£^2=  Kir,  (D-11) 

where  K  =  0,  ±1,  ±2,  ....  Rearranging  this  equation, 
the  location  of  the  grating  lobes  is  described  by: 


t  -KV 
^p  AX 


(D-12) 


The  relative  power  at  the  peaks  can  be  found  by  substitut¬ 
ing  Eq.  (D-12)  into  (D-10): 

P  (f  )  =  M2A2  .  (D-13) 

a  ^p 


The  first  mul  of  Eq.  (D-10)  occurs  at: 


(Mff  AX£  )/X/  =  ?r  , 


(D  -14) 


Ul|  MM.w, 
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and  therefore  the  width  of  the  grating  lobes  (from  the  peak 
to  the  first  null)  is: 


XT 

MAX  * 


(D-15) 


There  are  M  -  1  equally  spaced  nulls  between  the  grating 
lobes,  occurring  whenever: 

(MffAX£)/V  =  mr  (D-16) 


where  i.  =  ±  1,  ±2, _ _  M-l;  therefore,  there  are  M  -  2 

sidelobes  between  the  grating  lobes. 

The  poaks  in  the  transform  pattern  can  be  steered 
along  the  £  ?xis  by  applying  a  linear  phase  slope  across 
the  array.  For  a  linear  phase  slope  of  p  cycles  across 
the  array,  the  phase  of  the  mth  element  is  defined  by: 


2nm 
M  P  ‘ 


(D- 17) 


Inserting  this  into  Eq.  (D-3),  the  transform  pattern  of  the 
array  is  now  described  by: 

M 

U  («.n)  =  A  I  exp(j  — E)  exp[-j|£(m  -^±)AX£]  .  (D-18) 
m  =  l  J 

Using  a  similar  analysis  as  given  above,  Eq.  (D-18) 
reduces  to: 


.  .  sin  [(MffAXQ/X/-  Tip] 

a'4'  sin  ((itAXO/X/  -  rrp/MJ  ' 

This  equation  has  a  peak  whenever: 


_  ,K  i  p  i  \r 

'  (K  M*  AX 


where  K 


0,  ±1,  ±2,  . .  . 


(D-19) 


(D-20) 
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This  pattern  has  the  same  shape  as  Eq.  (D-9)  ex¬ 
cept  that  now  the  peaks  are  shifted  along  the  £  axis  by 

P(MAic).  ^^erefore*  ^e  location  of  the  peaks  on  the  £ 

axis  is  proportional  to  the  phase  slope  across  the  array. 

Using  the  Rayleigh  resolution  criteria  (that  the 
peak  of  one  beam  be  in  the  null  of  an  adjacent  beam),  it 
can  be  shown  from  Eqs.  (D-15)  and  (D-12)  that  there  are 
M  resolvable  beam  positions  between  the  ambiguities. 

Whe  i  the  input  function  is  a  two-dimensional 
orthogonal  periodic  array  of  isotropic  radiators  (as  shown 
in  Fig.  10),  Eq.  (D-l)  can  be  written  in  the  series  form 
(ignoring  a  eonstand  phase  term): 


N 

M 

)-  £ 

£ 

n  =  l 

m=l 

where 

M 

N  ■ 

AX 

AY 

j 


(D-21) 


tion,  and 

the  spacing  of  the  radiators  in  the  Y  direc¬ 
tion. 

Using  a  method  similar  to  that  used  to  obtain  Eq. 
(D-9)  from  Eq.  (D-3),  the  pattern  in  the  transform  plane 
for  this  two-dimensional  array  is  found  to  be 


Ua(£,r?)  =  A 


sin  [(MffAX£)/X/I 
sin  [Iff  AX£)/X/i 


sin  [(NffAYT?)/X/] 
sin  [(itAYtj)/X^ 


(D-22) 


This  equation  describes  a  two-dimensional  or¬ 
thogonal  periodic  array  of  grating  lobes  whose  peaks 
occur  whenever: 
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,  _K  X/ 
4  ‘  AX 


LX/ 

V  =  "ay  * 


(D-23) 


where  K  =  0,  ±1,  ±2,  ...»  and 
L  =  0,  ±1,  ±2,  ...  . 

The  widths  of  these  grating  lobes  in  the  £  and  r\  direc¬ 
tions  are: 


=  and 

ATJ  r  NAY  * 


(D-24) 


X.J  u  l 

The  relative  power  of  the  peaks  is  proportional  to  M2N2A2, 

The  peaks  in  the  transform  pattern  can  be  steered 
in  two  dimensions  by  applying  a  two-dimensional  phase 
slope  across  the  array.  Consider  the  following  two- 
dimensional  phase  slope  where  the  phase  of  the  (mfn) 
element  is  described  by: 


„  _  /  2rrm\  /  2irn\ 

^mn  “\m/P  \  N  /  * 


(D-25) 


The  plane  of  constant  phase  is  tilted  such  that  there  are  p 
complete  cycles  of  phase  across  the  array  in  the  X  di¬ 
rection  and  q  complete  cycles  of  phase  across  the  array 
in  the  Y  direction.  Substituting  Eq.  (D-25)  into  Eq. 
(D-21),  and  assuming  A(m,  n)  =  A: 


N  M 


d,V)  =  A  Z  Z  exp[j  2 ,r<~£  +  ^>)]  exp[-i  +  nAYr,)] . 

n=l  m  =  l  ^ 


(D-26) 
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Again  using  the  same  analysis  as  before,  the  pattern  in 
the  transform  plane  is  described  by: 


U  <£,tj)  =  A 


sin  l(MffAX£)/A/  -  rrp] 

sin  [<Nff  AYt)>/A/  -  irq} 

sin  l<ffAX$)/A/-”frp7M] 

sin  [(ffAY^)/X/  -  ffq/N) 

(D-27) 


This  equation  is  similar  to  Eq.  (D-22)  except  that  the  peaks 
are  shifted  along  the  £  axis  by  p(A//MAX),  and  along  the 
7)  axis  by  q(A//NAY).  Thus  the  position  of  the  grating 
lobes  is  a  function  of  the  two-dimensional  phase  slope 
across  the  array. 

Again  using  the  Rayleigh  resolution  criteria,  Eqs. 
(D-23)  and  (D-24)  can  be  used  to  show  that  there  are  M 
resolvable  beam  positions  between  the  grating  lobes  in 
the  £  direction  and  N  resolvable  beam  positions  in  the 
7}  direction  for  a  total  of  M  x  N  grating  lobes  in  the  two- 
dimensional  area  between  the  grating  lobes. 

So  far  the  discussion  has  been  limited  to  arrays  of 
isotropic  radiators.  In  actuality,  most  radiators  are 
directional  and  therefore  the  equations  described  so  far 
must  be  modified  by  the  "antenna  pattern"  of  the  individual 
radiators.  When  the  radiators  have  identical  transform 
patterns,  the  transform  pattern  for  the  array  is  de¬ 
scribed  by  the  equation: 

U(£,t?)  -  U  (£,rj)  U  (£,tj)  ,  (D-28) 

€  cl 

where  U(£,tj)  =  the  transform  pattern  for  an  array  of 
nonisotropic  radiators, 

U  (£,tj)  =  the  transform  pattern  of  the  array  when 
isotropic  radiators  are  assumed,  and 

U  (£,tj)  *  the  transform  pattern  of  the  individual 
radiator  elements. 
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When  0  volts  is  applied  to  the  masked  Membrane 
Light  Modulator  elements,  the  elements  appear  as  uni¬ 
formly  illuminated  circular  apertures.  As  shown  in  Ref. 
10,  the  transform  pattern  of  a  uniformly  illuminated 
circular  aperture  is  the  familiar  Airy  disk  which  is  de¬ 
scribed  by  the  Bessel  function: 


Ue(P,<P> 


Jr  ’ 

2?rpR> 

'  X/ 


(D-29) 


where  R  =  the  radius  of  the  mask  apertures  (20  jim),  and 

^  2*  ...2.  2.J 

p  +v  f  ■ 

The  first  zero  of  this  element  pattern  occurs  at 
2rroR  =  1.22*  ,  or 

(D-30) 


X/ 


0  = 


_  0.61  \f 


R 


The  transform  pattern  of  the  MLM  array  is  found 
by  substituting  Eqs.  (D-29)  and  (D-22)  into  Eq.  (D-28) 
(remembering  that  for  the  MLM,  AX  =  AY  and  M  =  N): 


U(C,T?)  *  C 


1  v  } 


sin  [(MffAXn)/X/l 
sin 


(D-31) 


From  Eq.  (D-30)  and  (D-23),  it  can  be  shown  that  along 
the  £  axis  there  are  seven  grating  lobes  under  the  main 
lobe  of  the  element  pattern. 
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When  voltages  are  applied  to  the  modulator  elements, 
they  deflect  as  paraboloids.  The  circular  mask  apertures 
are  now  considered  to  be  illuminated  by  a  parabolic  wave- 
front.  The  phase  across  the  elements  is  circularly  sym¬ 
metric  and  is  described  by  (Eq.  1,  Section  2): 

<p<r>=<po[l  -<|f>2]  . 

o 


where  *  the  peak  phase  at  r  =  0,  and 

R  =  the  radius  of  the  modulator  elements, 
o 

The  transform  pattern  for  the  masked  modulator  elements 
is  now  described  by 


R  2tr  /  r  2  \ 

=  J  J  exp(-jpo[l  -  (— )  ]  j  exp( - j  —  pr  sin  8)  r  dr  d0  .  (D-32) 

O  O  .  >  O  / 


R  2 

ue(fl>  -  i  «*P<  -j<PQ)J  exP[+jP0  >  1  J0(^)  r  dr  .  (D-33) 


R  2a 

'  J 

o  o  f 

Equation  (D-32)  can  be  reduced  to; 

R  2 

r 

o  o 


This  equation  has  not  been  solved  in  closed  form. 
However,  the  theoretical  variation  of  intensity  at  bore- 
sight  (p  =  0)  as  a  function  of  the  peak  phase  can  be  de¬ 
termined.  The  light  amplitude  on  boresight  is  given  by; 


Ue(0) 


P.  r  2 

e*P<-jP0)J  exp  [+j<P0(p->  1  r  dr 


expiW0  (^) 


(D-34) 
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From  Eq.  (D-34),  the  power  at  boresight  can  be 
shown  to  be: 

P(0)  *  U(0)u*(0) 


[l-coso  (A  ]  (D-35) 

O  rv 

O 

Equation  (D-35)  is  plotted  in  Fig.  18  for  RQ  =  50  and 
R  =  20  |im. 


=  C 


ttR' 


o 


<P 


o 


The  multi-channel  crystal  phase  modulator  can  be 
considered  as  a  linear  array  of  rectangular  apertures. 

As  shown  in  Ref.  10,  the  transform  pattern  of  a  uniformly 
illuminated  rectangular  aperture  is  described  by: 


ueU.r?)  =  ab 


,ira£  “ 
s.n  (^1 

'sinO~ 

(££l) 

V  J 

J 

(D-36) 


where  a  =  the  width  of  the  rectangular  aperture  parallel 
to  the  X  axis,  and 


b  =  the  width  parallel  to  the  Y  axis. 


Substituting  Eqs.  (D-36)  and  (D-9)  into  Eq.  (D-28\ 
the  transform  pattern  for  the  crystal  is  described  by: 


U(£,77;  *  Aab 


Bin  (if a. j/\f) 
(ira£/X/) 


sin  (irbrj/Xf)  a  in  (MttAX^/X/) 
(ffbrj/i/)  sin  (ffAX£/X/) 


(D-37) 


For  *he  present  crystal  mask,  =  b  =  AX/5.  From  Eq. 
(D-37)  it  can  be  shown  that  along  the  £  axis  there  are  nine 
grating  lobes  under  the  main  lobe  of  the  element  pattern. 
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APPENDIX  E 


THE  EFFECT  OF  A  SMALL  MAGNIFICATION  CHANGE 
ON  THE  FOCUS  POSITION  FOR  A  LENS  NEAR  UNITY 
MAGNIFICATION 


It  will  be  shown  in  this  appendix  that  the  total  ob¬ 
ject  plus  image  distance  for  a  lens  can  be  considered  con¬ 
stant  when  the  lens  is  used  near  unity  magnification.  The 
procedure  required  for  aligning  the  transmit  mask  and 
imaging  lens  is  greatly  simplified  by  this  fact  because  once 
the  image  is  approximately  at  unity  magnification,  the 
transmit  mask  no  longer  has  to  be  moved  (along  the  optical 
axis)  and  the  exact  magnification  can  be  obtained  by  moving 
just  the  imaging  lens. 

The  simple  lens  equation  is: 


l 

/ 


p  q 


where  p  =  object  distance. 


q  =  image  distance,  and 


(E-l) 


f  =  focal  length. 

If  we  let  the  object  plus  image  distance  equal  S, 
and  the  magnification  equal  m: 

S  =  p  +  q  ,  (E-2) 

m  =  q/p  .  (E-3) 

Inserting  Eqs.  (E-2)  and  (E-3)  into  Eq.  (E-l)  and 
rearranging  yields: 


o  =  f(m+l)2 
*  m 


(E-4) 
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The  minimum  of  this  function  can  be;determined 
by  setting  dS/dm  =  0  and  is  found  to  occur  at  m  =  1.  Since 
the  slope  is  zero  at  this  point,  the  value  of  S  around  this 
point  is  fairly  constant  for  small  changes  of  m.  There¬ 
fore  small  changes  in  the  magnification  around  unity  mag¬ 
nification  can  be  made  without  a  change  in  the  iinage-plus- 
object  distance. 

I 

i 

The  range  of  magnification  over  which  this  approxi¬ 
mation  is  valid  for  the  MLM  system  is  calculated  as  fol¬ 
lows.  Assume  that  a  change  in  the  magnification  of  Am 
results  in  a  change  of  the  object  plus  image  distance  of 
AS  (and  vice  versa).  Equation  (E-4)  becomes: 


2 

%/(m+Am  +  2  +  m-Am  +  Am  ) 
ai  /(2m+2  +  Ani2)  . 


At  m  =  1,  S  =  4/,  so  that: 

4/ +  AS«3  /(4  +  Am2)  ,  or 
2 

AS  as  /Am  * 


(E-5) 


The  depth  of  focus  for  a  diffraction  limited  lens 
in  a  coherent  system  when,  the  image  is  in  the  focal  plane, 
is  given  as  (Ref.  10): 


AS  as  ±  2 


(f)\t 

\d/  A* 


For  the  MLM,  the  image  is  at  2/  and  this  changes  to 


AS*,  ±  8 


(O’- 


(E-6) 
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Since  the  f  number  of  the  imaging  lens  is  // 4.  8,  and  the 
illuminating  wavelength  is  X  =  6328  k,  Eq.  (E-6)  becomes: 

AS  es  0.  12  mm  . 

Therefore,  the  percentage  change  in  magnification 
that  can  be  made  without  going  out  of  focus  is  (from  Eq. 
(E-5)  with /  =  150  mm): 


.  2’  0.12  mm  0  ,  A-4 

Am  =  -r-rrr -  =  8x10 

150  mm 


or 


Am  %  ±  0.  03 
Am  «  3%  . 


In  practice,  it  was  found  that  when  the  image  was 
within  approximately  1%  of  the  correct  magnification,  the 
transmit  mask  position  remained  constant. 
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APPENDIX  F 

SELECTION  OF  THE  ELEMENT  LOCATIONS  FOR  THE 
THINNED  MLM  ARRAYS 


Since  the  elements  along  each  row  of  the  MLM 
array  are  controlled  by  a  single  electrode,  the  phase 
of  each  of  these  elements  is  the  same.  Phase  slopes 
parallel  to  the  electrode  structure  are  required  in  order 
to  perform  two-dimensional  beam  steering,  and  there¬ 
fore  only  one  element  from  each  row  can  be  used.  A 
mask  is  used  to  block  the  light  from  all  but  the  appro¬ 
priate  element. 

The  criteria  to  be  used  for  selecting  the  location 
of  the  elements  in  the  thinned  array  are: 

1.  The  resolution  in  the  transform  pattern  of 
the  entire  100  x  100  array  should  b\*  main¬ 
tained. 

2.  The  sidelobe  levels  in  the  transform  pattern 
should  be  as  low  as  possible. 

If  independent  access  could  be  made  to  each  element  in 
the  100  x  100  array,  the  beam  could  be  steered  into 
10  000  resolvable  positions  (in  a  100  x  100  array)  using 
the  Rayleigh  resolution  criteria  which  requires  that  the 
peak  of  one  beam  occur  in  the  null  of  an  adjacent  beam. 
When  the  locations  of  the  thinned  array  elements  are 
chosen  judiciously,  this  resolution  of  the  whole  array  is 
maintained.  When  only  100  elements  are  used,  there 
will  not  necessarily  be  distinct  nulls  associated  with  the 
main  lobe  and  therefore  some  other  resolution  criteria 
is  required.  The  intensity  at  the  crossing  point  for  two 
adjacent  Rayleigh  beam  positions  is  approximately  4  dB 
down  from  the  intensity  at  the  peak,  and  therefore,  a 
4-dB  beamwidth  will  be  used  as  the  "Rayleigh"  resolution 
criteria  for  the  thinned  array. 


Preceding  page  blank 
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In  order  to  simplify  the  following  discussion,  the 
nomenclature  shown  in  Fig.  10  will  be  used.  For  the  MLM 
array,  Y  corresponds  to  the  direction  parallel  to  the  elec¬ 
trode  structure  (and  therefore  X  is  perpendicular  to  elec¬ 
trodes),  Also,  M  =  N  =  100  and  AX  -  AY  =  50  ^m.  In  the 
transform  plane,  £  is  the  direction  for  one-dimensional 
beam  steering  using  the  whole  100  x  100  MLM  array.  For 
simplicity,  isotropic  radiators  will  be  assumed. 

In  order  to  avoid  the  observation  of  false  beam  posi¬ 
tions,  it  is  necessary  to  minimize  the  peak  power  in  lobes 
other  than  the  main  grating  lobes.  The  minimum  sidelobe 
level  that  can  be  achieved  is  limited  by  the  requirement 
that  only  one  element  occur  on  each  electrode.  The  trans¬ 
form  pattern  along  the  £  axis  is  described  by  a  (sin  £)/£  dis¬ 
tribution  in  the  vicinity  of  the  main  grating  lobe.  The  peak 
intensity  of  the  first  sidelobe  for  this  distribution  is  13  dB 
down  from  the  peak  intensity  of  the  main  lobe.  The  goal 
therefore  is  to  keep  the  peak  intensity  at  any  other  point 
in  the  area  between  the  grating  lobes  below  this  -13  dB 
level. 


So  far,  no  restriction  has  been  made  on  the  loca¬ 
tion  of  the  elements  in  the  Y  direction.  Let's  consider 
what  would  happen  if  the  Y  location  for  each  row  were 
the  same.  This  results  in  just  a  linear  array  in  the  X 
direction.  The  resolution  of  the  100  x  100  array  is  main¬ 
tained  in  the-  £  direction  but  there  is  no  resolution  in  the 
rj  direction.  In  order  to  maintain  the  resolution  of  the 
whole  array,  the  width  of  the  lobe  along  the  rj  axis 
should  be  the  same  as  along  the  £  axis.  If  the  require¬ 
ment  is  imposed  that  every  column  contain  one  element 
(as  well  as  every  row),  then  the  transform  pattern  (in  the 
vicinity  of  the  main  lobe)  along  the  rj  axis  is  a  (sin  ?7)/r; 
distribution  (same  as  along  the  £  axis)  and  therefore 
the  width  of  the  main  grating  lobe  along  the  £  and  rj 
axes  satisfy  the  resolution  criteria.  This  does  not  mean 
however  that  the  resolution  of  the  array  is  necessarily 
maintained  in  all  directions.  A  good  example  of  this  is 


-  126  - 


TM  JOHN#  HOMtlNft  UNiVOKITV 

APPLIED  PHYSICS  LABORATORY 
S<vw«*  mm.  harwjmo 


the  case  where  the  100  elements  are  chosen  along  the 
diagonal  of  the  100  x  100  array  (which  again  reduces  to  a 
linear  array).  The  width  of  the  lobe  along  the  £  and  rj 
axes  is  correct,  but  there  is  no  resolution  along  the  £  =  rj 
line  and  therefore  there  are  only  100  resolvable  positions 
(in  the  direction  of  the  line  £  =  -r j). 

Probably  the  simplest  arrangement  of  the  thinned 
array  elements  would  be  a  set  of  straight  lines.  The 
transform  pattern  of  the  thinned  array  can  be  determined 
by  multiplying  the  transform  patterns  of  these  straight 
lines  (or  linear  arrays).  It  can  be  shown  that  for  linear 
slopes  at  the  input,  the  power  between  the  grating  lobes 
is  concentrated  in  line  patterns  (the  line  patterns  of  the 
linear  arrays).  Thus,  it  is  difficult  to  obtain  low  side- 
lobes  using  straight  line  patterns.  In  order  to  obtain  low 
peak  sidelobes,  it  is  necessary  to  spread  this  power  as 
uniformly  as  possible  over  the  areas  between  the  grating 
lobes. 


A  "random”  orientation  of  the  elements  is  re¬ 
quired  in  order  to  obtain  this  "thumbtack"  response.  In 
choosing  a  pseudo-random  orientation  for  the  100  elements, 
the  elements  are  still  restricted  to  one  element  for  each 
row  and  column  in  order  to  maintain  the  resolution  along 
the  £  and  rj  axes. 

A  computer  program  has  been  written  that  calcu¬ 
lates  the  power  at  any  point  in  the  transform  pattern  for  a 
thinned  two-dimensional  array.  Since  the  power  in  the 
transform  plane  is  symmetrical  about  the  origin,  the 
power  is  computed  over  only  half  the  unambiguous  trans¬ 
form  area.  The  number  of  pcints  at  which  a  calculation 
should  be  made  is  a  function  of  the  highest  spatial  fre¬ 
quency  component  (or  smallest  sidelobe  width)  that  is  ex¬ 
pected  in  the  transform  pattern.  It  can  be  shown  that  a 
total  of  80  000  points  satisfy  the  Nyquist  sampling  criteria. 

The  computer  program  itself  can  be  used  to  per¬ 
form  some  of  the  analysis  of  the  transform  pattern. 
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Threshold  levels  can  be  set  on  the  output  data.  A  plot  of 
the  power  in  the  transform  plane  can  be  obtained.  (A  one¬ 
digit  code  corresponding  to  the  relative  power  is  printed 
at  each  sample  point  which  crosses  the  threshold  level.  ) 
The  (£,77)  coordinates  of  the  sample  points  which  cross 
the  threshold  and  the  relative  power  at  the  points  can  also 
be  printed  on  a  separate  sheet. 

The  amount  of  computer  time  required  to  calculate 
the  power  at  the  80  000  samole  points  (using  the  IBM  360/ 
91)  is  approximately  3  minutes  (a  cost  of  around  $50). 

The  threshold  level  determines  the  number  of  points  to  be 
plotted.  An  additional  2  minutes  of  computer  time  and 
1  hour  of  plotting  time  is  required  for  a  typical  plot  with 
the  threshold  set  at  -15  dB,  resulting  in  a  total  cost  of 
around  $90.  Although  the  computer  plot  represented  an 
ideal  method  for  displaying  the  data,  only  a  few  were 
made  because  of  the  high  cost.  Since  the  peak  sidelobes 
were  typically  around  -10  dB,  the  threshold  was  set  at 
around  -11  dB  and  the  (£,77)  coordinates  and  power  of  the 
points  which  crossed  the  threshold  were  printed  out. 

These  points  (typically  10  to  20  points)  were  then  plotted 
by  hand. 

Several  pseudo-random  input  patterns  have  teen 
generated.  For  each  of  these  input  patterns,  the  power 
in  the  transform  plane  was  calculated.  The  best  of  these 
patterns  still  had  several  sidelobes  which  were  only  10 
dB  down  from  the  main  grating  lobes. 

A  limited  effort  was  made  to  "optimize'1  the  loca¬ 
tions  in  the  input  plane  in  order  to  eliminate  the  high  lobes 
in  the  transform  patterns.  A  high  sidelobe  in  the  2-D 
transform  plane  is  an  indication  that  a  particular  spatial 
frequency  in  a  given  direction  has  a  predominate  compo¬ 
nent  in  the  input  pattern.  The  particular  frequency  and 
direction  can  be  determined  from  the  (£>*7)  coordinates 
of  the  lobe.  When  the  input  plane  was  inspected,  a  pre¬ 
dominance  of  this  frequency  was  observed.  The  locations 
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of  the  elements  in  the  input  pattern  were  adjusted  in  order 
to  reduce  the  predominance  of  this  frequency.  Typically, 
this  resulted  in  the  position  change  of  around  10  elements. 
The  transform  pattern  for  this  new  input  pattern  was  then 
computed.  As  expected,  the  original  lobe  was  always 
eliminated.  However,  each  time  this  was  tried,  several 
new  sidelobes  would  pop  up.  This  "optimization"  proce¬ 
dure  was  attempted  on  three  different  patterns  and  each 
time  the  original  pattern  had  fewer  high  sidelobes. 
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